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CHAPTER 1 


INTRODUCTION 


1.1 Motivation 

Planetary gears are widely used in the transmissions of helicopters, automobiles, 
aircraft engines, heavy machinery and marine vehicles. Figure 1.1 illustrates a single- 
stage planetary gearset consisting of a sun gear, a ring gear, several planets, and a carrier. 
Any of the carrier, ring, and sun can be selected as the input or output component, and the 
power is transmitted through multiple paths of the planet meshes. Planetary gears have 
substantial advantages over parallel shaft drives, including compactness, large torque-to- 
weight ratio, diminished loads on shafts bearings, and reduced noise and vibration due to 
the relatively smaller and stiffer components. 

Despite planetary gears’ advantages, noise and vibration remain major concerns in 
their applications. In most helicopters, planetary gears are used in the last stage of gear 
reduction. This planetary gear is mounted directly to the helicopter cabin, so its vibration 
is the main source of structure-borne cabin noise, which can exceed 100 dB (Krantz, 
1993). Measurement of the cabin noise shows that gear mesh frequency and its harmonics 
are the dominant acoustic frequencies (Figure 1.2). Extensive cabin noise results in 
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Figure 1.1 Planetary gear diagram (Lynwander, 1983) 



Frequency, kHz 


Figure 1.2 Helicopter cabin noise spectra (Heath and Bossier, 1993) 
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operator fatigue, communication difficulty and health risk from extended exposure. 
Mitigating the cabin noise becomes a crucial requirement for effective operation of 
military and civilian helicopters. In vehicle automatic transmissions, planetary gear 
vibration can be transferred through bearings and mounts to the passenger compartment 
and results in discrete high frequency pitches in automotive interior noise. Consumers 
often perceive the gear noise as poor quality and mechanical problems. Vibration 
reduction in planetary gears provides substantial benefits; reduced noise, improved 
reliability, more efficient power transfer, and reduced maintenance costs. However, the 
noise and vibration issues in current planetary gear design rely mostly on empirical 
experience and cut-and-try methods, rather than engineering understanding. General 
design guidelines are needed to minimize planetary gear vibration. 

Vibration reduction of planetary gears requires thorough examination of the 
structural dynamics. A detailed review of planetary gear dynamic analyses before 1980s 
is found in August (1983). These analyses on planetary gear dynamics include modeling 
and free vibration investigations (Cunliffe et al., 1974; Botman, 1976; Frater et al., 1983; 
Velex and Flamand, 1996), neutralization of transmission error excitation (Seager, 1975), 
dynamic tooth load due to run-out errors (Hidaka, 1979e; 1980), mesh stiffness variation 
(Kasuba and August, 1984), load sharing among planets (Ma and Botman, 1984), and 
torsional vibration and dynamic loads (August and Kasuba, 1986). More analytical 
studies were performed after 1990s. Kahraman derived a nonlinear, time-varying planar 
model (1994a) and subsequently extended it to three-dimensions and examined the 
influence of planet phasing on dynamic response (1994b). Using a three-dimensional 
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model for helical gears, Kahraman and Blankenship (1994) investigated the load sharing 
and mesh phasing among planets. Kahraman (1994c) also reduced his model to a purely 
torsional one to predict natural frequencies and vibration modes. More recently, Agashe 
(1998) and Parker et al. (2000a) used a finite element tool to investigate the dynamic 
response and planet phasing issues in planetary gears. This special computational tool 
naturally includes the time-varying mesh stiffness and transmission error without ad hoc 
specification of these factors. Parker (2000) also rigorously proved the effectiveness of 
using planet phasing schemes to suppress planetary gear vibration. 

The experimental studies on practical planetary gear vibrations are scarce due to the 
difficult access to the internal gears. Chiang and Badgley ( 1 973) investigated the noise 
spectra generated from ring gear vibrations in the planetary reduction gearbox of two 
helicopters (Boeing-Vertol CH-47 and Bell UH-1D). Toda and Botman (1979) 
experimentally showed that planetary gear vibration resulting from spacing errors can be 
minimized by proper indexing of the planets. Botman (1980) presented some typical 
measurement results on the planetary gear of a PT6 aircraft engine. His experiments 
showed some peculiar behavior of planetary gear vibration regarding load sharing, 
response due to gear errors, and dynamic instability. Hidaka and his colleagues 
experimentally studied the dynamic behavior of a Stoeckicht planetary gear and 
published a series of reports (Hidaka et al, 1976a,b; 1977; 1979a~d). Their reports 
studied some important issues such as load distribution, effect of different meshing-phase 
among sun/ring-planet meshes, etc. Velex et al. (1994) matched the natural frequency 
measurement of a Stoeckicht epicyclic train with their finite element calculations. Rakhit 
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(1997) measured the subsynchronous vibrations at the turbine bearings in a gas 
turbogenerator and proposed a design of the epicyclic gearbox to reduce the vibration. 
Kahraman (1999) developed a generalized model to predict load sharing of planets under 
quasi-static conditions and validated the model with experiments. 

According to a comprehensive literature research, less analytical investigations 
have been done on planetary gear dynamics than those for parallel shaft gears. This is 
largely due to the modeling complexity of planetary gears. Important complications 
include multiple mesh contacts, detailed kinematics, mesh stiffiiess variation, 
transmission error excitation, contact loss nonlinearity, elastic ring gear vibration and 
geometric imperfections. Most previous research uses numerical or experimental methods 
to examine specific planetary gear systems. Some critical issues remain unsolved and 
require systematic analytical study. 


1.2 Critical Topics 

The fundamental task of analytical planetary gear research is to build a dynamic 
model. For different analysis purposes, there are several modeling choices such as a 
simple dynamic factor model, compliance tooth model, torsional model, and geared rotor 
dynamic model (Ozguven and Houser, 1988). According to the source-path-receiver 
relationship between the planetary gear, bearing/mounting, and the cabin, different 
boundaries can be selected for building the model. This study focuses on the 
understanding of planetary gear dynamic behavior, so a single stage gearset with discrete 
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elements is the basic model for investigation. In previous lumped-parameter models 
(Cunliffe et al., 1974; Botman 1976; Kahraman, 1994a. b, c), the gyroscopic effects 
caused by carrier rotation have not been considered. Because planetary gears have planets 
mounted on the rotating carrier, the Coriolis and centripetal accelerations caused by the 
carrier rotation introduce gyroscopic terms into the system model. For high-speed 
applications such as aircraft engines (30,000 rpm), gyroscopic effects may heavily impact 
the system stability and behavior. This project will develop a planetary gear model 
including the gyroscopic effects, contact loss nonlinearity, mesh stiffness variation, and 
static transmission error excitation. Despite the use of the term planetary gear, this model 
is applicable for epicyclic gears with any configuration (fixed/floating sun, ring, and 
carrier, and non-equally spaced planets). The model is the fundamental tool for the 
analytical research. 

The free vibration analysis calculates critical parameters such as natural frequencies 
and vibration modes that are essential for almost all dynamic investigations. Cunliffe et 
al. (1974) numerically identified the planetary gear natural frequencies and vibration 
modes for a specific thirteen degree-of-ffeedom system. Similar work has been done by 
Botman (1976), Frater et al. (1983), and Kahraman (1994c) for other example planetary 
gears, but no systematic characterization has been obtained. This project reveals the 
unique structure of natural frequency spectra and vibration modes due to the cyclic 
symmetry of planetary gears (Lin and Parker, 1999a). All the vibration modes are 
classified into rotational, translational and planet modes with distinctive properties. The 
structured vibration modes are rigorously characterized for general epicyclic gears and 
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validated by the computational results from a finite element model. These well-defined 
properties are not valid when planets are arbitrarily spaced, but still apply to practically 
important case of diametrically opposed planets (Lin and Parker, 2000a). The free 
vibration properties are very useful for further analyses of planetary gear dynamics, 
including eigensensitivity to design parameters, natural frequency veering, planet mesh 
phasing, and parametric instabilities from mesh stiffness variations. 

Another key issue is how design parameters affect the natural frequencies and 
vibration modes. During the design process, model parameters are often altered to 
evaluate alternative design choices, reduce weight, and tune the system frequencies to 
avoid resonance. The influence of design parameters on planetary gear natural 
frequencies was touched on in a few papers, but general conclusions were not presented. 
In the plots of natural frequencies versus design parameters, veering phenomena (Leissa, 
1974; Perkins and Mote, 1986) often occur and obstruct the tracing of eigenvalue loci 
under parameter changes. In the veering neighborhood, where two eigenvalue loci 
approach each other and then abruptly veer away, the veering vibration modes are 
strongly coupled and change dramatically (Figure 1.3). It is necessary to systematically 
study natural frequency and vibration mode sensitivities and their veering characters to 
identify the parameters critical to planetary gear vibration. In addition, practical planetary 
gears may be mistuned by mesh stiffness variation, manufacturing imperfections and 
assembling errors. For some symmetric structures, such as turbine blades, space 
antennae, and multi-span beams, small disorders may dramatically change the vibration 
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Figure 1 .4 Stability chart of Mathieu equation x + (S - 2f cos2 t)x = 0 (Meirovitch. 

1970). The hatched areas are instability regions. 

modes and result in mode localization (Pierre, 1988; Cornwell and Bendiksen, 1992; 
Happawana et al., 1998). Vibration modes with dominant motion localized in one planet 
lead to load sharing unbalance, which can severely undermine the power transfer 
efficiency and damage the gear teeth and bearings. This work presents a thorough 
eigensensitivity analysis of the natural frequencies and vibration modes to key model 
parameters for both tuned (cyclically symmetric) and mistuned planetary gears. The 
parameters considered include mesh/support stiffnesses, component masses, moments of 
inertia, and operating speed. Taking advantage of the structured vibration mode 
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properties, natural frequency sensitivities are expressed in simple, closed-form formulae 
which relate the sensitivities to modal strain and kinetic energy (Lin and Parker, 1999b). 
Well-defined veering rules are derived from these formulae and vibration mode 
properties to predict veering and its strength. The influence of design parameters is 
examined through a benchmark example. The knowledge of natural frequency spectra, 
vibration mode properties, eigensensitivity formulae, and special veering rules are 
combined to provide considerable insight into planetary gear free vibration. 

It is well known that mesh stiffness variation is a major excitation source of gear 
vibration. For spur gears, the time-varying stiffness is caused by the alternating number 
of teeth in contact. It is a periodic function at mesh frequency, which is the number of 
tooth mesh cycles per second. The mesh stiffness variation serves as a parametric 
excitation and results in instability under certain conditions. Parametric instabilities are 
particularly dangerous because they can occur at excitation frequencies well below the 
system natural frequencies. The mesh stiffness variation can be further complicated by 
the interaction of transmission error excitation (Smith, 1983) and contact loss 
nonlinearity (Blankenship and Kahraman, 1995; Kahraman and Blankenship, 1996, 
1997). Literature reviews of parametrically excited systems can be found in the work of 
Ibrahim and Barr (1978). For a single pair of gears excited by harmonically varying 
stiffness, Bollinger and Harker (1967) used the one degree-of-freedom Mathieu equation 
to determine the instability regions. The instability conditions are often illustrated in the 
plots of the exciting frequency versus the amplitude of varying stiffness, as shown in 
Figure 1.4. Benton and Seireg (1978, 1980a) obtained the response to mesh stiffness 
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variation and external excitations at integer multiples of the rotation speed. They 
experimentally verify the resonance region obtained from simulation and demonstrated 
the damage of parametric instability on gear teeth. Amabili and Rivola (1997) studied the 
steady state response and stability of the single degree of freedom system with time- 
varying stiffness and damping. Other researchers (Benton and Seireg, 1981; Kahraman 
and Blankenship, 1996; Nataraj and Whitman, 1997; Nataraj and Arakere, 1999) also 
investigated gear parametric instabilities using single degree-of- freedom models. For 
multi-mesh gear systems, it is surprising to find little investigation on parametric 
instability in the published literature. Although Benton and Seireg (1980b, 1981) studied 
a gear system with two meshes, they uncoupled the model into two single degree-of- 
freedom equations with some simplifications. Their conclusions on the mesh stiffness 
phasing effect contradict another investigation (Tordion and Gauvin, 1977) using an 
infinite determinant analysis (Bolotin, 1964). This conflict will be clarified using 
perturbation analysis (Hsu, 1963, 1965; Nayfeh and Mook, 1979) and numerical 
integration methods. In addition, this work extends parametric analysis in two-stage gear 
systems to planetary gears where parametric excitations are more complicated as 
different contact ratios and phasing conditions exist between the sun-planet and ring- 
planet meshes. August and Kasuba (1986) and Velex and Flamand (1996) numerically 
computed dynamic responses to mesh stiffness variations for planetary gears with three 
sequential phased planets. Their results showed the dramatic impacts of mesh stiffness 
variation on dynamic response, tooth loading, and load sharing among planets. The 
operating conditions leading to planetary gear parametric instability have not been 
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analytically investigated. The well-defined vibration mode properties will be used to 
derive the operating conditions leading to planetary gear instability. The boundaries are 
expressed in simple forms and the effects of contact ratios and mesh phasing are 
analytically investigated. In practical design, planet mesh-phasing schemes are often 
applied to cancel or neutralize the excitations from transmission errors (Seager. 1975; 
Kahraman, 1994a; Kahraman and Blankenship, 1994; Parker, 2000). This study shows 
that particular instabilities are eliminated under certain phasing conditions, which can be 
achieved by selection of design parameters according to the analytical results. Dynamic 
response and tooth separations induced by parametric instability are numerically 
examined. 


1.3 Scope of Investigation 

The scope of this project is to advance the modeling and understanding of planetary 
gear dynamics and analytically examine certain critical factors affecting planetary gear 
noise and vibration. This research focuses on the analytical investigation of the following 
specific tasks. 

• Derive a lumped-parameter model for spur planetary gears, including different planet 
phasing, gyroscopic effects, mesh stiffness variation, and transmission error 
excitation. The model is valid for general epicyclic gears with any number of planets 
and will be the fundamental tool for further research. 
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• Analytically characterize the unique structure of the natural frequency spectra and 
vibration modes of general planetary gears. The cases with equally and arbitrarily 
spaced planets (including diametrically opposed planets) are considered. 

• Use the vibration mode properties to obtain simple, closed-form formulae for the 
eigensensitivities to important design parameters. According to these formulae and 
characterized natural frequency veering rules, the effects of design parameter changes 
on planetary gear free vibration are investigated. 

• Investigate the parametric instabilities caused by multiple time-varying mesh 
stiffnesses. Two-stage gear systems are examined first to clarify pervious conflicts and 
derive simple expressions of instability boundaries. Then, the analytical method is 
extended to planetary gear systems. The well-defined modal properties are used to 
identify the effects of contact ratios and mesh phasing on planetary gear parametric 
instability. 
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CHAPTER 2 


PLANETARY GEAR MODELING 


2.1 Modeling Considerations 

Lumped-parameter modeling is used in this project for dynamic analysis. All gears 
are considered as rigid bodies and component supports are modeled by springs. A single- 
mesh model is shown in Figure 2.1 for a pair of spur gears. All bearings/supports are 
modeled as two perpendicular springs. The gear tooth meshes are represented by springs 
acting on the line of action with parallel viscous dampers. The transmission error e(t), a 
prescribed displacement input, is included as indicated. The tooth separation nonlinearity 
acts as a step function h(t) where h=J when the teeth are in contact and h=0 when the 
teeth lose contact. This model can be extended to planetary gears with multiple meshes. 

In planetary gear modeling, the following factors must be considered. 

1 . Assumptions. The analysis deals with planar vibration of single stage planetary gears. 
Helical gears require three-dimensional modeling and are not considered here. 
Excellent lubrication is assumed and tooth friction forces are neglected. No damping 
is included, although a viscous damper could easily be added in parallel with the 
mesh and bearing springs. 
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Figure 2. 1 A single mesh model 

2. Versatile configurations. There are many configurations for general epicyclic gears. 
By fixing one or more of the coaxial components (carrier, ring, and sun), various 
configurations such as planetary, star, and differential gears can be obtained. In 
typical designs, one of the coaxial members is free to translate to enhance the load 
sharing among planets. The number of planets and their positions also vary in 
practical applications. The proposed model should accommodate these configurations 
and be flexible for general applications. 

3. System coordinates. Several choices of coordinates may be used in the modeling. The 
frame can be fixed or rotating; the planet coordinates can be parallel to each other or 
use local radial and tangential directions. The coordinate selection does not change 
the physical properties of the system, but may greatly affect the analysis difficulty. 


N AS A/CR— 200 1-210939 


16 



2.2 System Equations 

The planetary gear model used in our analysis is shown in Figure 2.2. Each 
component has three degrees of freedom: two translations and one rotation. The model is 
similar to that used by Kahraman (1994a) with two distinctions: (1) the planet deflections 
are described by radial and tangential coordinates, and (2) gyroscopic effects induced by 
carrier rotation are modeled. The radial and tangential coordinates more naturally 
describe the vibration modes. Gyroscopic effects in high-speed applications such as 
aircraft engines may dramatically alter the dynamic behavior from that at lower speeds. 
The coordinates illustrated in Figure 2.2 are used. The carrier, ring and sun translations 
Xh, }’h, h=c,r,s and planet translations £ n ,t] n ,n = l,-- N are measured with respect to a 
rotating frame of r reference i, j, k fixed to the carrier with origin o. The x h , h=c,r,s are 
directed towards the equilibrium position of planet 1, and C, n , rj n are the radial and 
tangential deflections of the n-th planet. The basis i, j, k rotates with the constant carrier 
angular speed £2 C . The rotational coordinates are u h = r h 8 h ,h = N , where 9h is 

the component rotation; r h is the base circle radius for the sun, ring and planet, and the 
radius of the circle passing through the planet centers for the carrier. Circumferential 
planet locations are specified by the fixed angles y/„, where y/„ is measured relative to the 
rotating basis vector i so that y/j = 0. The details of the model derivative are given in Lin 
and Parker (1999a). The equations of motion are 

Mq + D f Gq + [K b + K„, - Q^K n ]q = T(r) + F (t) (2.1) 

q = (s c , y c ,u c ,x r ,y r ,u r ,x s ,y s , u s , ^, ih,u l; - 

carrier ring sun planet 1 planet N 
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where the matrix components are given in the Appendix A. M is the inertia matrix and 
K b is the bearing stiffness matrix. G and Ko result from high-speed carrier rotation and 
have not been included in published models. To model the time-varying stiffness 
associated with changing numbers of teeth in contact at each mesh, K m can be 
decomposed into mean and time-varying components. Tooth separation nonlinearity is 
implicitly included in K m (t). T (t) is the applied external torque and F (t) represents the 
static transmission error excitation. 

This comprehensive model is the fundamental tool for further research of planetary 
gear dynamics. It is readily applicable for specific configurations of epicyclic gears. 
Assigning a large value to a transverse or torsional stiffness restrains the translation or 
rotation of a component; assigning a very compliant stiffness in a component support 
floats the component. 



Figure 2.2 Planetary gear model 
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CHAPTER 3 


NATURAL FREQUENCY AND VIBRATION MODE PROPERTIES 


This chapter analytically investigates the natural frequency spectra and vibration 
modes of general planetary gears. The cyclic symmetry of planetary gears leads to 
highly-structured free vibration characteristics that are identified herein. Unique 
properties of the eigensolutions for the linear time-invariant case are presented for an 
example system. The identified properties are then mathematically shown to be 
characteristics of general planetary gears. 

3.1 Equally Spaced Planets 

The free vibration of the linear, time-invariant representation is 

Mq + 0 c Gq + (K-fi c 2 K Q )q = O (3.1) 

where K=K b +K m . The carrier speed is assumed to be small and the gyroscopic terms G 
and Kq are neglected. The associated eigenvalue problem of (3.1) is 

Q)f = K(pi (3.2) 

where a>i are natural frequencies and <j> : = [p f? p r ,p ;r ->Pi> , ”PAr] r are vibration modes with 
Pz, = [x h ,y h ,n h ] T ,h = c,r,s for the carrier, ring and sun, and p„ = [£„,%,u„] T for the 
planets. At this stage, the planets are assumed identical and equally spaced; all planet 
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bearing stiffnesses are equal k p „=k p , all sun-planet mesh stiffnesses are equal (k x „=k sp ); 
and all ring-planet mesh stiffnesses are equal km-krp. With these specifications, planetary 
gears are cyclically symmetric structures that can be divided into N identical sectors. The 
cyclic symmetry of planetary gears leads to distinctive natural frequency and vibration 
mode properties that will be demonstrated analytically. 

We first illustrate the eigensolution properties through a numerical example with 
the parameters in Appendix B case I. Typical vibration modes for N=4 are shown in 
Figure 3.1, where the movements of the carrier and ring are not shown in order to clarify 
the figures. Some interesting conclusions are obtained from Figure 3. 1 : 

1. Six natural frequencies always have multiplicity m=l for different N. Except for the 
zero natural frequency, their values increase as additional planets are introduced. 
Their associated vibration modes have pure rotation of the carrier, ring and sun 
(Figure 3. la), so these modes are named rotational modes. 

2. Six natural frequencies always have multiplicity m=2 for various N. Some natural 
frequencies increase monotonically while others decrease monotonically as N 
increases. The carrier, ring and sun have pure translation in the corresponding 
vibration modes (Figure 3. lb), so these modes are defined as translational modes. 

3. Three natural frequencies have multiplicity m=N-3 and exist only if N>3. Their 
associated vibration modes are termed planet modes because the carrier, ring and sun 
do not move; only planet motion occurs in these modes (Figure 3.1c, d, e). For each 
of these three natural frequencies, the corresponding vibration modes span an N-3 
dimensional eigenspace. 
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This example shows there are at most fifteen different natural frequencies for N>3; 
additional planets only change the multiplicity of the planet mode natural frequencies. 
When N<3 , all natural frequencies are distinct and the vibration modes do not have 
special structure because of the loss of cyclic symmetry. In nearly all planetary gear 
designs, three or more planets are used to take advantage of the load sharing among 
planets and subsequent discussion is restricted to this case. The eigensolution properties 
identified in the example are analytically shown to be true for general planetary gears 
(Lin and Parker, 1 999a) and summarized below. 

Planet Mode. A planet mode has the form 

<t>, = [0,0,0, Wjp,,— w*p,] r (3.3) 

where w„ are scalars (vr/=/) satisfying 

X sin y/„ =0 X u '« cos ¥„=° X u '» =0 ( 3 - 4 ) 

Planet modes have the following characteristics: 

(i) The associated natural frequency has multiplicity N-3, 

(ii) The translation and rotation of the carrier, ring and sun are zero. 

(iii) The planet deflections are scalar multiples of the first (or any other arbitrary) 
planet's deflection components. 

Rotational Mode. A rotational mode has the form 

<P, =[Pr.Pr'Ps.Pl>- "Plf 

Rotational modes have the following characteristics: 

(i) The associated natural frequency is distinct. 
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(ii) The translation of the carrier, ring and sun are zero, i.e., p h = [0,0, u h ] 7 , where j=c,r,s , 


(iii) All planets have identical deflections, i.e., p, =p 2 = - =p v = [£■,./;, .uj 7 
Translational Mode. A pair of translational modes has the form 


<P,= 


Pc 


Pc 

p, 


P, 

p s 

(cos^,p, -l-sin^p,) 

<P,= 

Pv 

(—sin ^,p, + cos ^,p, ) 

(cos^ A ,p, +sin^p,) 


(-sin^ A .p, +cos^p,) 


(3.6) 


Translational modes have the following characteristics: 

(i) The associated natural frequency has multiplicity two. 

(ii) The rotation of the carrier, ring and sun are zero. Furthermore, the carrier, ring, and 
sun translations in the degenerate modes $ and $ are related by p h = [x h ,y h ,0i\ T and 

Pa = [-v„,*a,0] 7 \ h=c,r,s, 

(iii) The planet deflections for a pair of vibration modes are related by 

(3-7) 

where p„ and p„ = are the deflections of the n-th planet in <p, and <£>, . I is 

a 3x3 identity matrix and y/ n = 2n(n - 1) / N . 


P„ 


cos y/ n I sin y/ n I 

Pi 

_P»_ 


-siny/„I cos^ 7 I 

Pi. 
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(a) rotational mode 


(b) translational mode 



(e) rotational planet mode 


Figure 3.1 Typical vibration modes. Dashed lines are the equilibrium positions 
and solid lines are the deflected positions. Dots represent the component centers. 
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3.2 Arbitrarily Spaced Planets 

In general, much of the above well-defined structure of the natural frequency 
spectra and vibration modes is lost when the planets are arbitrarily spaced. A notable 
exception is the planet modes. Additionally, for the practically important case of 
diametrically opposed planets, the free vibration retains its unique properties. 

Planet Mode: As for planetary gears with equally spaced planets, systems with 
arbitrary planet spacing always have three sets of planet modes of the form (3.3) with 
multiplicity N-3. Only the coefficients w„ obtained from (3.4) are affected by y/ n . 

Rotational Mode: In general, the rotational and translational modes couple together 
for arbitrary planet spacing and no special modal structure can be identified. For certain 
planet spacing, however, they still have distinguishing properties. A case of particular 
interest is that of diametrically opposed planets, which is common in industrial 
applications. Consider a system with each of N/2 pairs of planets located along arbitrarily 
oriented diameters. A pair of opposing planets have the position relation y/ n +N 2 = Wr> +71 hi 
this case, 

Z sin yr„ = 0 Z co Wn =0 Z w » =0 (3.8) 

Accordingly, systems satisfying (3.8) have six rotational modes with property (3.5). For 
arbitrarily distributed planets not satisfying (3.8), rotational modes do not exist. 

Translational Mode: While the translational modes couple with the rotational 
modes for truly arbitrary planet spacing, they retain their structure for systems satisfying 
(3.8). The notable difference with equally spaced planet systems is that the natural 
frequencies are no longer degenerate because the cyclic symmetry is disturbed. To start 
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with, the planet deflection relations in a translational mode are derived from (3.7) for any 
three planets i,j, k, 

sin(^, ~Wj)Pk +sin(^ ; -^)p, +sin(^ t -^)p, = 0 (3.9) 

Thus, the n-th planet deflection can be expressed as a linear combination of pi and p 2 - 
The component modal deflections for a translational mode become 

P/, =[x h ,y h ,u h ], h = c,r,s, p„ =[sin(^ 2 -^„)p 1 +sin^,p 2 ]/sin^ (3.10) 

Therefore, planetary gears with planet positions satisfying (3.8) (for example, 
diametrically opposed planets) have twelve distinct vibration modes that have the special 
structure (3. 10) of a translational mode. 


3.3 Modal Strain and Kinetic Energy 

The vibration modes can be further characterized by modal strain and kinetic 
energy. The total modal strain energy U and kinetic energy T are related to the natural 
frequencies and vibration modes by 

U = \ <t>J K4 =U C + U cu +U r + U„ +U S + U su +f J (U„+U ni + U sn ) (3.11) 

2. n=\ 

T =1+ T„ + l + T m +T, + T m +f,(T,+T.„) (3.12) 

2 M=1 

where (i) Uh, Uhu, h= c,r,s are the strain energies in the translational and rotational 
support springs, respectively, of the carrier, ring and sun; (ii) U„, U nh U sn , n=l, .... N are 
the strain energies in the n-th planet bearing, ring-planet mesh and sun-planet mesh; and 
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(iii) Th , Thu • h=c, r, s and T„, T mi are the modal translational and rotational kinetic 
energies, respectively, of the carrier, ring, sun, and planets. The detail of these individual 
energies are 

= 2 K (*?, + yl \ U /,„ = ~ K, u\, h = c , r, s 


ey(x; + v,;). 

nr hi _ 2.2 

hi ~ ^ 2 fy 11 h !* 

h = c,r,s 

(3.13) 

7 2 / r2 2 \ 

-a, (£ +n„x 

nn ^ P 2 2 

T »u = - 2 «„ 

2r p 




where the mesh and bearing deflections 6 are defined below. 

sun-planet mesh: S m =y, cosy/ sn -x, sin^ OT -£ sina t - 77 ,, cosor, + + «„ + e„, (3.14) 

ring-planet mesh: 5 m =y r cos^ m -x r sin^„, sin« r + 7 „ cosa r + «, -w„ +e m (3.15) 

planet bearing radial: £ m . = y r sin y/„ + x. cos^„ - £ n (3.16) 

planet bearing tangential : S nl = y c cos y/ n - x c sin (3.1 7) 

For the example planetary gear (Appendix B, case II), the modal strain and kinetic energy 
distributions are shown in Figure 3.2, and the dominant motion and strain energy are 
listed in Table 3.1. According to (3.13), the dominant kinetic energy occurs in the 
component with dominant motion. Modal strain energy indicates those vibration modes 
most susceptible to parameter variation and identifies the most heavily loaded component 
for response. In the next chapter, the modal energy distribution will be used in the 
derivation of simple formulae to calculate eigensensitivity to key design parameters. 
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3.4 Discussions 


Planet modes of multiplicity N-3 are remarkably insensitive to planet location and 
retain their special properties for arbitrary planet spacing. Coupling between rotational 
and translational modes occurs for arbitrary planet spacing, and distinct properties can 
not be identified. For systems satisfying (3.8), however, rotational and translational 
modes have structured properties. This includes the common case of equally spaced and 
diametrically opposed planet pairs. The foregoing development applies to general 
epicyclic gear configurations. Configurations having fixed or floating components are 
obtained by letting the associated bearing stiffness approach zero or infinity to obtain the 
eigensolution properties. For example, for the cases of a fixed ring, fixed sun or fixed 
carrier, a 3(N+2) degree of freedom system is obtained. The vibration modes in such 
systems consist of five rotational, five pairs of translational and three groups of planet 
modes. 

The free vibration properties have been validated through a finite element 
computation. A finite element model (Figure 3.3) for this system was built by Agashe 
(1998) and Parker, et at. (2000) and analyzed with CALYX (Vijayakar, 1991). 
Computational modal analyses were performed by applying impulse inputs and obtaining 
the frequency response functions. The identified natural frequencies match the analytical 
results (Table 3.2) within 3.3 percent difference. 
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(a) Normalized modal strain energy at the sun/planet meshes (1-4), ring/planet meshes 
(5-8), planet bearings (9-12), carrier and sun bearings (13, 14). 
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(b) Normalized modal kinetic energy at the carrier, sun, planet translations (1,2, 3-6), and 
carrier, sun, planet rotations (7,8,9-12). 


Figure 3.2 Modal energy distributions of the system in Appendix B, case D. 
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No. 

/(Hz) 

Type 

Dominant motion 

Dominant strain energy 

1 

0 

R 

Rigid body mode 

2.3 

825 

T 

All components 

Carrier, sun bearings 

4 

1661 

R 

All components 

Planet bearings 

5 

1808 

P 

Planet radial translation 

Planet bearings 

6,7 

1834 

T 

Sun, planet translation 

Sun, planet bearings 

8 

1985 

R 

Sun rotation, planet translation 

Planet bearings 

9,10 

2326 

T 

Sun, planet translation 

Sun, planet bearings 

11 

5964 

P 

Planet tangential translation 

All meshes 

12,13 

6429 

T 

Planets 

Ring-planet meshes 

14 

6451 

R 

Planets 

Ring-planet meshes 

15 

6982 

P 

Planet rotation 

All meshes 

16,17 

10430 

T 

Sun translation 

Sun-planet meshes 

18 

13068 

R 

Sun rotation 

Sun-planet meshes 


Table 3.1 Dominant motion and modal strain energy in the vibration modes of the 
system in Appendix B. R: rotational mode, T: translational mode, P: planet mode. 


Mode Type 

Trans 

Rot 

Plan 

Trans 

Rot 

Trans 

FEM (Hz) 

778 

1144 

1729 

1676 

1723 

2110 

Analytical (Hz) 

769 

1156 

1609 

1710 

1781 

2175 

Difference (%) 

1.2 





-3.0 


Table 3.2 Comparison of analytical and FEM natural frequency analyses. 
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Figure 3.3 Finite element model of the planetary gear 


The well-defined eigensolution properties are useful for subsequent research 
because almost all analytical investigation of planetary gear vibration phenomena 
ultimately require the natural frequencies and vibration modes. This includes, for 
example, critical dynamic behaviors such as forced response to static transmission error, 
use of planet phasing to eliminate excitation of particular modes, parametric instability 
from time-varying mesh stiffness, contact loss nonlinearities, natural frequency and 
vibration mode sensitivity to key design parameters and gyroscopic effects. The specific 
properties characterized in this paper theoretically explain the selective participation of 
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the different classes of vibration modes in planetary gear dynamic response (Parker et al., 
2000a; Parker, 2000). Additionally, the structured vibration mode properties lead to 
simple, exact formulae to calculate natural frequency and vibration mode sensitivity to 
parameter changes (Lin and Parker, 1999b). Using these properties, well-defined veering 
rules of natural frequency loci are analytically derived (Lin and Parker, 2000b). These 
special veering rules help to trace the evolution of the loci for changes in the design 
parameters and identify the veering zones where vibration modes undergo dramatic 
changes. Apart from analytical applications, the identified natural frequency structure 
provides important information for tuning the system frequencies to avoid resonance. 
There are at most 15 different natural frequencies in general planetary gears with N 
planets, so there are only 15 potential resonant frequencies. 
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CHAPTER 4 


EIGENSENSITIVITY TO DESIGN PARAMETERS 


During the design process, system parameters are varied to evaluate alternative 
design choices, avoid resonances, optimize load distribution, and reduce weight. It is 
important to characterize the effects of parameter variations on the natural frequencies 
and vibration modes for effective vibration tuning. In planetary gear dynamic models 
(Figure 2.2), the key design parameters include the mesh stiffnesses, support/bearing 
stiffnesses, component masses, and moments of inertia. The influence of some design 
parameters on planetary gear natural frequencies was touched on in a few papers. Botman 
(1976) and Cunliffe et al. (1974) both presented plots of natural frequencies versus planet 
bearing stiffness. Kahraman (1994c) showed the effects of mesh/bearing stiffnesses on 
the natural frequencies in his torsional model of planetary gears. Saada and Velex (1995) 
discussed the influence of ring support stiffness on free vibration. These analyses were 
based on parametric studies of example planetary gears and assume the planetary gears to 
be cyclically symmetric (tuned) systems. Eigensensitivity analysis for mistuned systems 
is necessary to identify the critical modes that are susceptible to irregularity. Frater et al. 
(1983) studied the vibration modes with one unbalanced mesh stiffness, but general 
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conclusions were not obtained. This project analytically investigates the natural 
frequency and vibration mode sensitivity to most system parameters such as 
mesh/bearing stiffnesses, component masses and moments of inertia. Simple, closed-form 
expressions are obtained to calculate eigensensitivities for both tuned and mistuned 
system, in addition, eigenvalue veering phenomena are investigated to identify dramatic 
changes of natural frequencies and strong coupling of vibration modes. Design guidelines 
are summarized from eigensensitivity and veering analyses to predict influences of 
system parameters on planetary gear free vibration. 


4.1 Calculation of Eigensensitivity 

The eigensensitivity analysis calculates natural frequency and vibration mode 
derivatives with respect to stiffnesses, masses, moments of inertia and the carrier rotation 
speed Q c . Eigensensitivity to stiffness and inertia design parameters are examined in the 
absence of gyroscopic effects Q c . Gyroscopic effects are important in high-speed 
applications such as aircraft engines, and eigensensitivity with respect to Q c is studied 
separately in section 4.4. The eigenvalue problem for the study is in form (3.2), i.e., 

(K-A,M)<*,=0 (4.1) 

where X,=ox. The eigensensitivity for problems in the form (4.1) has been thoroughly 
investigated (Courant and Hilbert, 1953; Adelman and Haftka, 1986; Friswell, 1996) and 
the necessary results are introduced below. The unique modal properties of planetary 
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gears are then invoked to reduce these general results to simple expressions specific to 
planetary gears. 

Let ( )' and ( )" denote the first and second derivatives with respect to a model 
parameter (i.e. mesh/bearing stiffness, component mass or moment of inertia). For 
simplicity, the eigenvalue derivatives and A," are calculated; the relations 6^'= A, 7(2(0,) 
and (Ot "= (2XjAj "-A, 2 )/ (4 co, 3 ) yield the natural frequency sensitivities. For a distinct 
eigenvalue, the eigensensitivities are (Fox and Kapoor, 1968; Rogers, 1970) 


4'=$ r (K'-4M')$ (4-2) 

1 _ , ^<P k T ( K'-A,M')A 

a>;= - T (0, r M + X , -i> k (4.3) 

^ *=1 A; A, k 

A,"= A,M')</>, + <t>! { K' - A.M' - 2k‘ ')</>, (4.4) 


For the case of degenerate eigenvalues, consider a system having a group of 
eigenvalues X with multiplicity m. The first-order eigenvalue derivatives A ,' are 
the eigenvalues of 

Da, = A'a t , D= ^(K'-^M'F (4.5) 


where T=[yi y m ] is an arbitrary set of independent eigenvectors associated with this 

degenerate eigenvalue and is normalized such that F r MT = I mxm . For the case when all 
Aj' obtained from (4.5) are distinct, the a, are uniquely obtained with the normalization 
aj T aj =1. This procedure determines the set of independent eigenvectors $= Ta, that admit 
continuous change of the eigenvectors as the degenerate eigenvalues split into distinct 
ones when a parameter is varied. 

The eigenvector derivatives for distinct are expressed as (Friswell, 1996) 
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<=V, +Oc,., 


/= 1 , -,m 


(4.6) 


where 


_ + 

' ^ 2 ; Vk 

k=m + 1 A A- 


(4.7) 

x'-x; 

j*i 

(4.8) 

c„ =--<P t M'® 
2 


(4.9) 


The second derivatives are 

X'= 2<l>l{K’-X i M')\ l +<t>J(K”-X,M'-2X'M')<t> l , /= l,-,m (4.10) 

For the case when all X, 'obtained from (4.5) are degenerate, the a, are not unique 
and hence m are arbitrary in the eigenspace. The eigenvector derivatives can 

not be determined when these degenerate modes do not separate. However, X," can be 
obtained from the eigenvalues of (Friswell, 1996) 

E = 2<t> T (K'-A,M')V + <t> T (K'-X i M'-2X'M')<& (4.11) 

where V=[vi,. . .,v m ] is determined by (4.7). X, "are not affected by the selection of O. 

The foregoing development is used subsequently to derive general, closed-form 
eigensensitivity relations for Xj', 0'and X, "for planetary gears. These expressions yield 
eigensolution approximations according to 

X = X + ^—\ P = p 0 {p- P q), 0 = 0 + 
p vp 

where p represents any system parameter with nominal value po and multiple parameter 
perturbations are permitted. Eigensolutions X, 0 are for a nominal set of model 


3~l p = Po(t Po) 
p dp 


(4.12) 
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parameters referred to as the unperturbed system, and the derivatives are evaluated for 
this unperturbed system. Eigensolutions I, 0 are for the perturbed system with varied 
parameters. The unperturbed system is taken to be tuned in this study. Note that this does 
not meaningfully restrict the results because parameter variations leading to both tuned 
and mistuned perturbed systems are examined. 


4.2 Eigensensitivity to Mesh and Support Stiffnesses 

The stiffnesses under consideration (Figure 2.2) include mesh stiffnesses km k sn , 
transverse support stiffnesses k c , k r , k s , k„ , and rotational support stiffnesses kcw inn k.su- 
The natural frequency sensitivity to a certain stiffness is found to be uniquely associated 
with the modal strain energy occurring in that spring. To demonstrate the procedure, let 
the sun-planet mesh stiffness k s „ be the varied parameter. 


4.2.1 Tuned System 


Considered the case where all sun-planet mesh stiffnesses k s „=k sp are altered 
equally so the perturbed system remains tuned. For rotational modes, the 
eigensensitivities are obtained from (4.2)-(4.4) 




dk 


sp 


n = 1 


_ y 1 y 1 KfiL . 

dk k=x „ =1 \ ~ \ 


(4.13) 

(4.14) 


k*i 
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(4.15) 




dk 


sp 


™K-K 

k*i 


<ZWL) ! 


H=1 


where S' sn is the spring deformation of the sun-planet n mesh in mode fa given (3.14). 
The rotational mode property (3.5) dictates that all sun-planet mesh deformations are 


equal, i.e.. = S‘ A , so (4.13) becomes 


d\ 


=N(s;y 


(4.16) 


In (4.14), fa'\s expressed as a modal expansion of eigenvectors, and the contribution of 
each eigenvector is readily obtained from the coefficients of fa. When two eigenvalues A, 
and Ak are nearly equal, the influence of fa on fa ' is dominant because the denominator in 
its coefficient is small. In such cases, the second derivative fa" is also large, and the 
natural frequency changes rapidly with k sp . 

The translational mode eigenvalues A, , do not separate because the perturbed 

system remains tuned. Thus, the matrix D in (4.5) has degenerate eigenvalues A/', A2' 
Accordingly, the unperturbed eigenvectors fa h2 can not be uniquely determined from the 

procedure associated with (4.5), and O=[0;, fa] are an arbitrary pair of translational modes 
of the unperturbed system. From (4.5), A\, A 2 ' are the eigenvalues of 


D = <f> r K'<D = 


Z 


0O 2 

jPJL 


(C) 2 _ 


(4.17) 
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Use of the translational mode property (3.6) yields y L _,(<?i)‘ ~ )~ - 

y N S l S 2 = 0. Thus, the eigenvalues of D (/.<?., A,' X 2 ) are degenerate and have the 

1 sn Sf} 

form (4.13) for i=l,2. From (4.1 1) and the translational mode properties, A\ ”, X 2 ” are 


d ~ _ y -1 y -1 2(S stl S sll ) 


dk 


*=3 n=l 




ij = h2 


(4.18) 


Planet modes are also degenerate and the procedure is similar to that for 

translational modes. For planet modes <$=[0/ 0?], the elements of the matrix D are 

ZX = for ij=l m. Applying planet mode property (3.3) to calculate 6' sn 

results in £>//=...= D mm and D 0 =0, It follows that all V of a group of planet modes 
are equal and can also be expressed as (4.13) for In the same way, all A, , 

are equal and of the form (4.18). 

Equation (4.13), which is valid for all three types of vibration modes, can be related 
to the modal strain energy U in 0,. With the definition of strain energy U sn in (3.13), 
(4.13) becomes 


d\ 

<*s P 



dco t 


p n = 1 


2X 


( 4 . 19 ) 


Equation (4.19) allows one to obtain the natural frequency sensitivity to sun-planet mesh 
stiffness by inspection of the modal strain energy distribution. 

As an example, consider a planetary gear used in the transmission of a U.S. Army 
OH-58 helicopter. The nominal model parameters are listed in Appendix B, case II. The 
natural frequencies from (4.1) are shown in Figure 4.1a for a range k sp . The strain 
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energies of each spring are calculated according to (3.1 1) and their distribution in mode 
16 (a translational mode) is shown in Figure 4.1b.c for two cases: k sp =70 Nl\xtn and 
k sp =500 Ni\im. The associated vibration mode 0k, is also shown for these two cases. 
Little strain energy is stored in the sun-planet meshes U sn for case 1, while substantial 
strain energy results in case II. Consequently, a>i6 is more sensitive to k sp in case II than 
in case I. This conclusion is consistent with the larger slope of the a>i 6 locus for case II in 
Figure 4.1a. In fact, natural frequency sensitivities to all stiffnesses can be obtained 
quantitatively directly from the strain energy distribution using (4.19) and analogous 
relations (4.20) and (4.21) below. 


d\ , 2 2 d\ 

= u l h = c,r,s 

K h 

(4.20) 

n= 1 K p w=l 

A 2 f 

njL 

UK rp n = 1 K rp n=\ 

(4.21) 


These relations apply for all three types of vibration modes. Expressions for <p , ' and A," 
for all of the stiffness parameters are collected in Lin and Parker (1999b). Recalling the 
special properties of vibration modes, (4.19)-(4.21) imply that 

1 . Rotational modes are independent of the transverse support stiffnesses of the carrier, 
ring, and sun because these components have no deformation of, and hence no modal 
strain energy in, their transverse support springs. Thus, use of a "floating" sun, ring, 
or carrier i.e., low stiffness support) has no impact on rotational modes. 

2. Translational modes are similarly independent of the rotational support stiffnesses of 
the carrier, ring, and sun. 

3. Planet modes are insensitive to all carrier, ring, and sun support stiffnesses. 
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Figure 4.1 (a) natural frequency versus the sun-planet mesh stiffness k sp . (b), (c) Mode 16 
strain energy distribution in case I and case II. All U are defined in (3. 11). 
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4.2.2 Mistimed System 


In practical planetary gears, mistuning may be caused by differing mesh stiffnesses 
between planets due to differing numbers of teeth in contact, manufacturing variations, 
and assembly errors. To study the effects of mistuning on eigensolutions, we examine the 
sensitivity to parameter variations that differ between the planets. Consider an example 
with only the first sun-planet mesh stiffness k s j varying from the nominal (unperturbed) 
value k sp . The derivatives of the mass and stiffness matrices with respect to k s j are 
M ,= M"=K"=0, K ' = oK / dk sX . 


The eigensensitivities of the rotational modes are obtained from (4.2)-(4.4) 


= (^) 2 = 


d\ 

* A 

dk _Zj 


r u * 

All 


si k 
k*t 


^ = t 


2 


dk 


<t> k 


(4.22) 

(4.23) 

(4.24) 


si 


KK-K 

k*i 

Equation (4.22) relates A/ to the modal strain energy in the first sun-planet mesh. 
Equations (4.22)-(4.24) are similar to (4.13)-(4.15) without the summation over n 
because the varying parameter is located only at the first sun-planet mesh. 


For translational modes, the eigensensitivities are (Lin and Parker, 1999b) 


\2 

— VVsl) ’ 


*,1 ™ ’ * MX 


= («) 2 =o 


d<t\ V * 


dksl 

K 


- = o 


0 


(4.25) 

(4.26) 

(4.27) 
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The behavior of fa and fa is shown in Figure 4.2 for the example of Appendix B. case II. 
A pair of translational mode natural frequencies separate as a disorder e=k s i/k sp -l is 
introduced. The modal strain energy distributions in the four sun-planet meshes are 
shown for e= 0, -0.1. fa is sensitive to k sI because of the high strain energy in the first 
sun-planet mesh, fa has no strain energy in the first sun-planet mesh and is independent 
of k sl . The linear ( A, + ek sp K ) and quadratic ( A, + £k sp % + £ 2 k sp K' ) approximations of 

the loci are shown in Figure 4.2 and agree well with the exact loci. These two loci 
intersect exactly at e = 0 when there is only one disorder in the perturbed system. If one 
more disorder £ 2 =k^k sp -l= 0.1 is added at the second mesh (Figure 4.3), the two loci 
suddenly change direction and veer away. For an initially tuned (cyclically symmetric) 
system, two independent varying parameters ( e.g k s i and k s j) are necessary to break the 
symmetry of both fa and fa and cause frequency loci veering (Happawana et al., 1998). 
The regular perturbation does not give a good approximation in the veering zone where 
two loci are close to each other. Singular or improved perturbation methods are needed 
(Pierre and Murthy, 1992; Wu, 1993; Chen et al., 1995; Lin and Lim, 1997). If the 
number of planets N= 4 or 5, the planet modes have multiplicity m=l,2 and their 
eigensensitivities are obtained from (4.22)-(4.24) or (4.25)-(4.27). When N>5, 
eigensolutions of matrix D in (4.5) are difficult to achieve in closed-form, but can be 
obtained numerically. 
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6.1 h 


6 1 1 1 1 1 1 5 ‘ i f 

0.2 0.15 0.1 0.05 0 0.05 0.1 0.15 0.2 0.25 0.3 

e= ksi/ksp -1 


Figure 4.2 Influence of the disorder e on the natural frequencies. Linear (...) and 
quadratic ( — ) approximations agree well with the exact loci ( — ■). 


x 10 4 



Figure 4.3 Influence of two disorders on the natural frequencies. Another disorder 
£ 2 =k s Vk sp -l =0. 1 at the second mesh is added to the system shown in Figure 4.2. 
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4.3 Eigensensitivity to Gear Mass and Inertia 


The parameters of interest consist of masses ( w c , m r , m s , ni p ) and moments of inertia 
(4,4, 4 I p ) for the carrier, ring, sun, and planets. When the perturbed system remains 
tuned, the eigenvalue derivatives for the three types of modes are 


n = 1 rn p n=l 


2 

- - — T 

l h 

1 N 

-~rL u n 
r p „=i 


(4.28) 

(4.29) 


where T h , T hu , h=c, r, s and T„, T nu are the modal kinetic energies defined in (3.12). 
Expressions for <p' and /l," for all of the mass and inertia parameters are collected in Lin 


and Parker (1999b). Figure 4.4a shows an example plot of the natural frequencies versus 
the sun moment of inertia I s . Most natural frequencies are insensitive to changes in I s . The 
kinetic energy distribution and vibration modes of mode 18 (a rotational mode) are 
shown in Figure 4.4b,c for cases I and n. The sun has more rotational kinetic energy T su 
in case I than in case II, so C0is locus has larger slope in case I. Equations (4.28) and 


(4.29) allow quantitative calculation of natural frequency sensitivity to all masses and 
moments of inertia directly from the modal kinetic energy distributions. Considering the 
properties of vibration modes, some conclusions are immediate from (4.28) and (4.29): 

1. Rotational modes are independent of the masses of the carrier, ring, and sun because 


such modes have no translations of these components. 

2. Translational modes are similarly independent of the moments of inertia of the 


carrier, ring, and sun. 

3. Planet modes are independent of both masses and inertias of the carrier, ring, and sun. 
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Tc Ts 


Tcu Tsu 



Tc Ts 
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Figure 4.4 (a) Natural frequency versus the sun moment of inertia /*. (b), (c) Mode 1 8 
kinetic energy distribution in case I and D. The T are defined in (3.12). 
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4.4 Eigensensitivity to Operating Speed 

In high-speed applications (e.g. aircraft engines), gyroscopic effects may 
significantly alter the system stability and dynamic behavior. Eigenvalue derivatives 
evaluated for Q c =0 are calculated to assess the influence of operating speed on the 
natural frequency spectrum. The gyroscopic eigenvalue problem of (3.1) is obtained from 
the separable solution q = 0 l e ja> ‘ 1 

[-a;, 2 M + jo),Q c G + (K - Q] K a M = 0 (4.30) 

For practical operating speeds £2 C =0 (i.e., subcritical), the eigenvalues remain purely 
imaginary. Suppose a zero speed natural frequency 0} has multiplicity m and the 
arbitrarily chosen independent eigenvectors are T=[j5, ..., y m ] with normalization 
r 7 Mr = l mx „, . While eigenvectors for £2 ^ 0 are complex, the y are real. Differentiation of 
(4.30) with respect to Q c and evaluation at £2 C =0 yield 

(K - co;M)<p' = (2a) i co'M - ja> i G)Tn i = f (4.3 1 ) 

where 0, =Ta, . Applying solvability and normalization conditions results in an mxm 
Hermitian eigenvalue problem 

Da, =oi, 'a,, D = yT r Gr/2 (4.32) 

The natural frequency sensitivities (o' are obtained from the eigenvalues of (4.32) for the 
three classes of vibration modes. Rotational mode natural frequencies are distinct and 
(4.32) becomes a scalar equation. Hence oi,'= 0 because yjGy t = 0 for real y, and skew- 
symmetric G. For translational modes yi and y 2 , D and its eigenvalues are 
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Figure 4.5 Natural frequency versus the carrier rotation speed 

/[ o y[Gy->l T 

D= J - * co' u =±ylGyJ2 (4.33) 

2[- y;Gy 2 0 

For a group of planet modes yi,...,y m , the properties (3.3) guarantee Z) ; = y/Gy ; = 0 for 

D„=0, i=l m because of the skew-symmetry of G. Thus, D=0 and all planet 

mode natural frequency sensitivities vanish, i.e., 0)'= 0. 

Equation (4.33) can be used to approximate the frequency loci a ), = co, + Q c 0)' . 
The result co'= 0 for rotational and planet modes at Q = 0 indicates the natural 
frequencies of these modes are scarcely affected by operating speed. Figure 4.5 shows the 
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first 10 frequency loci versus Q c for the gear system in Appendix B case II. The 
rotational mode (ft!/, ftfc) and planet mode ((ft) loci are nearly flat lines and not sensitive 
to operating speed. Translational mode frequencies (tftj, (ft, 7 , (ft,/o) split as Q c is 
increased from zero. In this example, (ft and (ft at Q c =600 rad/s deviate about 10 percent 
from the zero speed value. Typical helicopter carrier speeds are less than 100 rad/s. For 
applications with high speed (e.g. turbofan and turboprop engine systems), heavy 
component masses, and compliant stiffnesses, the gyroscopic effects can be more 
significant. If a natural frequency locus has large slope and decreases to zero in the range 
of operating speed, the stability and system behavior are dramatically impacted. 

4.5 Natural Frequency Veering 

In the plots of natural frequencies versus design parameters, eigenvalue veering 
occurs, where two eigenvalue loci approach each other as a parameter is varied but then 
abruptly veer away like two similar charges repelling (point B in Figure 4.6a). The 
phenomenon has been studied extensively (Leissa, 1974; Perkins and Mote, 1986; Pierre, 
1988; Chen and Ginsberg, 1992). The vibration modes of the veering eigenvalues are 
strongly coupled and undergo dramatic changes in the veering neighborhood. In the case 
of especially sharp veering, it is sometimes difficult to distinguish between intersection 
and veering just by observing eigenvalue plots. When multiple curves veer or intersect 
close together (Figure 4.6a), strong modal coupling, and the associated operating 
condition response changes that occur, are not identifiable from frequency loci plots. The 
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objective of this work is to analytically characterize the rules of eigenvalue veering in 
planetary gear vibration. Simple rules emerge to predict if two eigenvalues veer or cross. 
The veering sharpness is also calculated. 


4.5.1 Veering/Crossing Criterion 


A method for detecting eigenvalue veering/crossing in general dynamic systems is 
developed by Perkins and Mote (1986). When two eigenvalue loci veer away, their loci 
curvatures indicate the abruptness of curve direction changes. Perkins and Mote 
estimated the loci curvature in the veering neighborhood using coupling factor 

2[*, r (K'- A,M')0 r ] 2 


= 2[0! (K'- x r M'xps 
Xr X -X, 


Xs = 


X, - X r 


(4.34) 


The coupling factors Xr- Xs approximate the local curvatures and the coupling strength. 
Figure 4.6b shows X~Xu ~ ~X >8 versus the varying parameter k sp for the veering loci (O 14 
and cois in Figure 4.6a. Notice the sharply changing vibration modes indicated in Figure 
4.6b. The two veering loci exchange mode shapes from point A to C, even though the loci 
do not intersect. The modes are strongly coupled at B and do not look like either of the 
veering modes just outside the veering zone. When a parameter is adjusted in the veering 
zone, the drastic changes in the vibration modes can greatly impact the operating 
condition dynamic response, tooth loads, load sharing, and bearing forces and possibly 
lead to mode localization. The degree to which individual modes are excited by dynamic 
mesh forces (/.<?., the modal forces) also changes dramatically as veering alters the 
modes. If the coupling factors are all zero, A,- and X s loci cross; otherwise, veering occurs. 
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x 10 



k sp (N/|im) 


Figure 4.6 (a) Natural frequencies versus the sun-planet mesh stiffness. Natural 
frequencies are numbered under the nominal conditions (dashed line) in Table C.l. Three 
pairs of veering are loci 14 and 18, (12,13) and (16,17), and 11 and 15. (b) Coupling 
factor x of 0014 and co/s- These two loci exchange mode shapes from point A to C. 
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4.5.2 Veering Patterns In Planetary Gears 

When applied to planetary (or any epicyclic) gears, the coupling factors reduce to 
particularly simple forms because of the unique structure of the vibration modes. In the 
frequency plots of tuned planetary gears, loci of the same type (rotational, translation, and 
planet modes) never cross each other but veer away when they come close. The modes of 
different types can cross each other but never switch modes through veering. This special 
veering pattern is analytically proved using the well-defined vibration mode properties 
(Lin and Parker, 2000b). The general veering/crossing patterns are summarized in Table 
4. 1 . This pattern is generally valid when stiffness or inertia parameters are varied. 

The special veering patterns are helpful to trace the evolution of eigenvalue loci and 
identify the effects of design parameters on planetary gear vibration. The planetary gear 
in a helicopter powertrain is used as an example. The nominal model parameters are 
given in Table C.l, case II. Table 3.1 identifies the mode type and where the dominant 
strain energy is in each mode. The natural frequencies are numbered at the nominal 
conditions (indicated by the dashed lines in Figure 4.6-1 1). 
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P2 
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x/v* 
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! 

V 

P3 

X 

X 

X/V* 

V 

1 

! 


Table 4.1 Veering (V) and crossing (X) patterns in planetaiy gears. X/V means crossing 
only occurs when k sp = krp and Or=cy. 
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Figure 4.7 Natural frequencies versus the planet bearing stiffness k p . The planet mode 1 1 
(pure tangential type PI) crosses the planet mode 15 (no tangential motion type P2). The 
nominal conditions (dashed line) are listed in Table C.l, case II. 

Mesh stiffnesses k sp , krp (Figure 4.6and 4.8) have little influence on the low natural 
frequencies O)j~a>io ■ This is because these modes are governed by bearing stiffnesses 
(Table 3.1) that are much smaller than the mesh stiffnesses. Modes 15-18 have large 
strain energy in the sun-planet meshes and are affected by k sp (Figure 4.6); modes 1 1-14 
have substantial strain energy in the ring-planet meshes and are affected by k^ (Figure 
4.8). When k sp is reduced from the nominal value, the changing C 0 j 5 ~( 0 is approach the 
G)i]~(Oi 4 loci. Because loci of the same type can not intersect, veering occurs between 
rotational modes 18 and 14, translational mode pairs (16,17) and (12,13), and planet 
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modes 15 and 11. Below the veering zones ( k sp <100 N//jm ), modes 11-14 are very 
similar to modes 15-18 above the veering zones ( k sp >800 N/juw). In the same way, one 
can predict the trend of frequency loci as krp is increased (Figure 4.8). Using the derived 
veering patterns and modal properties, the actual modes affected by varying parameters 
can be detected easily although the plots are complicated by veering phenomena. 

Support stiffnesses kh, kh u , h=c,r,s of the carrier, ring, and sun can vary over a wide 
range depending on the configuration (fixing or floating these components). Rotational 
and planet modes are independent of the transverse support stiffness kh because they have 
no translation of the carrier, ring, and sun; only translational modes are affected by 
changes in kh (Figure 4.9a,b). Considering the veering effects, kh significantly affects only 
one pair of translational modes with dominant strain energy in the transverse supports. 
When the rotational support stiffnesses kh„ are altered, similar results are obtained (Figure 
4.9c,d), except it is the rotational modes that are susceptible to khu variations. 

Altering planet parameters affects most natural frequencies as all modes involve 
planet deflections, in general. Applying the derived veering results, five pairs of veering 
are identified for changing planet bearing stiffness in Figure 4.7: rotational modes 8 and 
18, 4 and 14, translational modes (9,10) and (16,17), (6,7) and (12,13), and planet modes 
5 and 15. For large stiffness k p >5000 N/jum, eight natural frequencies increase rapidly to 
outside the range of interest. Planet mass m p and moment of inertia I p also have 
significant influence on the natural frequencies (Figure 4.10). 
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4.6 Discussion and Summary 


Eigenvalue loci veering also occurs when degenerate modes of symmetric system 
are separated by small disorders. In planetary gears, the cyclic symmetry can be broken 
by differing mesh stiffnesses at each planet mesh, manufacturing variations, and 
assembly errors. For cyclically symmetric or periodic systems with small disorders and 
weak structural coupling, mode localization often accompanies eigenvalue loci veering 
(Pierre, 1988). Planetary gears have relatively strong coupling between the planets 
through the carrier and teeth meshes, so strong mode localization is unlikely even in the 
presence of loci veering. After examination of many cases with various configurations 
and parameters, the authors have not found a realistic example of mode localization in 
planetary gears. 

The special veering patterns of planetary gear eigenvalue loci are easily 
summarized. Two approaching eigenvalue loci of the same type (rotational mode, 
translational mode, and planet mode) veer away while two loci of different mode types 
cross each other. The mode shapes are exchanged across the veering zone. In the veering 
zone, the modes are strongly coupled and markedly different than outside the zone. One 
can expect significant differences in response from these changed modes. These rules 
result from planetary gears' unique modal properties and apply to all design parameters. 

The effects of key design parameters are summarized below: 
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1 . Mesh stiffness. k sp and k rp each control three different natural frequencies associated 
with one rotational mode, one pair of translational modes, and one group of planet 
modes. Dominant strain energy occurs in the tooth meshes of these vibration modes. 

2. Carrier, ring, and sun parameters, kh and h=c,r,s each affect only one natural 
frequency. The transverse stiffness kh controls one pairs of translational modes and the 
torsional stiffness kh„ controls one rotational mode. Floating or fixing the carrier, ring, or 
sun has limited influence on planetary gear modal properties. The carrier, ring, and sun 
masses and moments of inertia affect the same frequencies as their corresponding support 
stiffness, though the frequencies vary in the opposite direction. 



Figure 4.8 Natural frequencies versus the ring-planet mesh stiffness krp. Loci 12-14 will 
veer with loci 16-18 when krp is further increased. The nominal conditions (dashed line) 
are listed in Table C.l, case II. 
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3. Planet parameters. Planet bearing stiffness and planet inertia are the most influential 
parameters and affect most natural frequencies. A stiff planet bearing can be beneficial 
for resonance tuning because it substantially reduces the number of natural frequencies 
in the lower frequency range that is commonly of most interest. 



k $ (N/nm) 



k su (N/nm) 



k c (N/nm) 



k cu (N/jun) 


Figure 4.9 (a), (b) kh,h=s,c only affect translational modes, (c), (d) k hu , h=s,c only affect 
rotational modes. The nominal conditions (dashed line) are listed in Table C.l. 
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CHAPTER 5 


PARAMETRIC INSTABILITY FROM MESH STIFFNESS VARIATION 


In previous chapters, the mesh stiffnesses are considered constant. Current work 
investigates the impacts of mesh stiffness variation on planetary gear dynamics. The 
variation of mesh stiffness gives rise to a parametric excitation. The linear, time-varying 
system studied herein is governed by a set of coupled Hill’s equation 

Mq + [K 0 +K v (0]q = F(/) (5.1) 

where Ko is the mean stiffness matrix and K v (t) is the variational part of the stiffness 
matrix. F(t) is the external excitation. K v (t) is periodic with frequency £2 (mesh 
frequency). Mesh stiffness depends on element compliance, tooth error, profile 
modification, transmitted torque, and the position of contact (Kasuba and Evans, 1981). 
When K v (t) is simple harmonic in t, (5.1) becomes a set of coupled Mathieu equations. 

The parametric excitation from K v (t) causes instability and severe vibration when 
harmonics of the excitation frequencies are close to particular combinations of the natural 
frequencies. Three types of instability are of most interest: (1) primary instability £2=2 
(2) secondary instability £2=C0p, and (3) combination instability £2=0},+ 6%. The objective 
of this study is to systematically determine the operating conditions leading to parametric 
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instabilities in planetary gears. As discussed in Chapter 1, the existing analysis on 
parametric instability is scarce, inconsistent, and incomplete for multi-mesh gears 
systems. A simple two-stage gear chain is analyzed first, and then the analysis is 
extended to planetary gears with multiple meshes. Perturbation methods are used in this 
study to determine the boundaries separating the stable and unstable regions. Floquet 
theory and numerical integration are used to validate the analytical findings. The effects 
of contact ratios and mesh phasing on parametric instability are quantitatively identified. 
The interaction of instability, tooth separation nonlinearity, and dynamic response is 
discussed. 

5.1 Two Stage Gear Systems 

5.1.1 System Model 

Two-stage gear trains have three-gear and four-gear configurations (Figure 5.1). 
The gear bodies and intermediate shaft connecting gears 2 and 4 are assumed rigid. The 
tooth meshes are modeled as linear springs with stiffnesses ku, ku ■ The anchored shaft is 
flexible with torsional stiffness k/o. The gears have base radii /*,, i=l ,2,3,4. Only 
rotational vibrations 9i, 0 !?, 63 are considered. The equivalent masses are m/= Ij/rj, nij= 
h/rf , and m 2 = where /, are the moments of inertia of the gears and their 

connected shafts, v=r/r? for four-gear trains, and v=7 for three-gear chains. The 
equivalent stiffnesses are k 0 = k L0 /rf, k/= ku, and k 2 = v'ku . The shaft/gear rotations are 
measured by the base radius deflections xi= Ojr/, x 2 = $ 2 f 2 , and xj=^ryv. The system 
stability is governed by the free vibration equation 
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Figure 5.1 Two-stage gear systems with (a) four gears and (b) three gears. k Ll> k L2 
denote mesh stiffnesses and Z 2 , Z 4 are numbers of gear teeth, kio is the torsional 
stiffness of the anchored shaft. 


Mq + [K 0 +K v (O]q = 0 (5-2) 

where q = [x,, x 2 ,x 3 ] T . M=diag(m /, m 2 , m 3 ) is the inertia matrix. The stiffness matrix is 
represented by a mean value Ko and a variational part K r (t) as 



1 

+ 

O 

k ,i 

0 



K\ 

0 

K 0 = 

*,i 

0 

k g i + Ki 

kg2 

k ,2 
k *2 _ 

K ,.(0 = 

K 

0 

*vi + Ki 

Ki 

Ki 

K-2 . 


where k gi and k v ,(t) are the mean and time-varying components of the mesh stiffnesses, 
k/t) = k p + kjt) (Figure 5.2). The variational parts k v ,(t) are periodic at the mesh 

frequency Q, and expressed in Fourier series as 
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/ = 1,2 


(5.4) 


MO = sin sO,/ + b- l) coss£l,t), 

;=i 

where 2k a , is the peak-to-peak amplitude of k v , (Figure 5.2). The mesh frequencies Q\ and 
f2j are related by £2j=RQ, where 7?=ZyZ^ and Z?, Z 4 are the numbers of teeth on gears 2 
and 4 (Figure 5.1a). Note R=l, £2j=Qi for three-gear systems (Figure 5.1b). Mesh 
stiffness variation is obtained through measurement, calculation, or simple specification 
{e.g., sinusoidal or rectangular wave). For spur gears, rectangular waves are often used to 
approximate the mesh stiffness alternating between n and n+1 pairs of teeth in contact 
(Kahraman and Blankenship, 1999a). In this study, the k vi are specified as rectangular 
waves with variational amplitudes k ai , periods T, =£2/2n, contact ratios c„ and phasing 
angles p{Tj (Figure 5.2b). Thus, 

a) !) = — — sin[/7T(c- ; - 2/?,)]sin(/^t,), b)' ) = — —cos [ln{c < -2p : )]sin(/^c,) (5.5) 

In In 

in (5.4) for / = 1, 2, •••. Without loss of generality, one can specify pi=0, p?=h ( h is 
called mesh phasing). In practice, the first three or four Fourier terms reasonably 
approximate the mesh stiffness variation. 

For the time-invariant case, the eigenvalue problem associated with (5.1) is 

'I' 

K o 0, = a>~M0 i . The vibration modes <p, are normalized as <1> MO=I with (fh, <p 3 ]. 

Applying the modal transformation q = «J>u and using (5.4), equation (5.2) becomes 

+ , (£><'> sin /Q.t + ^cos/Q,/) 

r=l l=\ (5-6) 

+ 2e 2 (F^ sin lQ 2 t + Gj t ' r } cos lQ. 2 t)]u r =0, n = 1,2,3 
where e i =k cn /k gi , i-1, 2 and the matrices D, E, F, G are 
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Figure 5.2 Modeling of mesh stiffnesses k,{t) = k g , +M0- c> are contact ratios, k gi are 
average mesh stiffnesses, and p,T, are phasing angles. 


a[ a[ 0 

D (/) = ^,<P r a[ a[ 0 <D 

0 0 0 

a\ a[ 0 

F (,) = k g2 Q> T a\ a\ 0 O 

0 0 0 


'0 0 0 ~ 

E (,) = k gl <t> T 0 b[ b[ <t> 
0 b\ b[ 

0 0 0 ' 

G (,) = k g2 < D r 0 b' 2 b' 2 <S> 

0 b[ b‘ 2 


(5.7) 


Equation (5.6) is a set of coupled Hill’s equations subjected to multi-frequency 
parametric excitations from two gear meshes. 
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5.1.2 Conditions of Parametric Instability 

Parametric instability depends on the frequency, amplitude, and shape of the 
parametric excitations. In gear systems, these factors are directly associated with the 
operating speed and gear design parameters such as contact ratio, facewidth. diametral 
pitch, pressure angle, material properties, and so on. The corresponding model 
parameters are the stiffness variation amplitudes k a i, k a 2 , mesh frequencies 42/, 42_->. 
contact ratios c/, c?, and mesh phasing h. 

The variation amplitudes k m are assumed small compared to the average mesh 
stiffnesses k gi , so £,= k a /k gl «J in (5.6). At this point, £■/=£?=£ is specified; the case of 
£/*£? is discussed later. Using the method of multiple scales, three different mesh 
conditions are examined. 

1. Three-gear Systems: Equal Mesh Stiffness Variations 

In three- gear systems (Figure 5.1b), the two meshes have the same mesh 
frequencies Q, = Q 2 = Q . We consider the case where the gear facewidths and material 
properties, which primarily determine mesh stiffness for a given number of teeth in 
contact, are such that the amplitudes of mesh stiffness variation are the same at the two 
meshes (£, = £ 2 =£). The contact ratios and mesh phasing are allowed to differ between 
the two meshes, however. In practice, the contact ratios are changed using center 
distance, diametral pitch, pressure angle, tooth addendum, and other parameters. The 
mesh phasing depends on the layout of the gears and the numbers of teeth. 
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The parametric instability when sQ is close to Cty+ait, is considered. Let 
IQ = 0) p + (O q + oe , where a is a detuning parameter. Using multiple scale method in 


(5.6), the boundaries of the instability regions are (Lin and Parker, 2000c) 

Q=i(», +a . f ±^) a«=-L-k < 5 - 8 > 




For single mode instabilities (p=q), equation (5.8) becomes 

n = ± y-HD% + C ) ! +<£"’+ G !? ) ! ■ w > 2 • • • 

l t p 


(5.9) 


As an example. Figure 5.3 shows the boundaries for primary (/=/), secondary 
(1=2), and combination (p#q, 1=1 ) instabilities. The parameters are given in Table 5.1 and 
c l= c 2 = 1.5, h- 0. Floquet theory, numerical integration to compute the fundamental 
matrix, and the fundamental matrix eigenvalues (Nayfeh and Mook, 1979) are used to 
determine the actual instability regions denoted by * in the figures. The first-order 
approximations from (5.8) agree well with the numerical solution, even when e is not 
small. In Figure 5.3, the instability region around Q=2c^ is much larger than that of the 
primary instabilities around 2coi and 2oh. This is explained by examining the vibration 

modes. From (5.8), the primary instability boundary slopes are governed by . 


Expansion of D, E, F, G in (5.7) yields the diagonal terms 


2„<0 


&£ = MA, + C = WA, + A,)^ 

C = MA, + A,) 2 A (0 . G p p = k g2 (<p 2p +<p 3p ) 2 b^ 


(5.10) 
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where <j>i p , <jh p , <p3 P are the gear rotations in mode Let 8i p =<pi p +<jh p =v/+x? be the 
relative deflection of the first mesh in mode Similarly, 82 P =<f>2 P +<p3 P =X2+xs represents 
the modal deflection in the second mesh. For the primary instability boundary around 
2cq„ insertion of (5.10) into (5.8) yields 


Ki = KMM 0 + KM? 


(5.11) 


The mesh deflections Si, & in each mode can be observed from the mode shapes (Figure 
5.4). The two meshes in (f>j are both in phase and have smaller Si, S2 than those of 0?, 
where the two meshes are both out of phase. Thus, A33 > A*, 1 * and the instability 

boundaries around 2 cos have larger slope than those around 2q)j. Mode <jh has one mesh in 
phase and the other out of phase, so the size of its primary instability regions is between 
that of 0/ and fy ( Aj] > A ( -j, > A 1 , 1 , 1 ). In addition, the mesh deflections in a vibration mode 


are related to the modal strain energy U 1 = k x S x /2, U 2 = k 2 S \ / 2 . Examination of (5. 1 1 ) 
shows that vibration modes with more strain energy in the meshes have larger instability 
regions and are more susceptible to parametric excitations. The above results apply for 
mesh stiffness variations of arbitrary shape. 

For mesh stiffnesses having rectangular waveforms, one can clearly identify the 
effects of contact ratios and mesh phasing on the instability regions. Use of (5.5) in (5.1 1) 
yields 

A PP = [*X sin 2 (c,/r) + k] 2 S 2p sin 2 (c 2 n) + 

~k gl k 2 Sl Sl sin(C];r)si n (c 2 7r) cos (c, -c 2 +2h)n ]{— — ) 2 ^ ^ 

■; -- - - /rry„ 
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Inertia 

ni]-l, m2~0.3 , m 3 =4.0 

Average mesh stiffness 

kgi~ kg?—! 

Shaft stiffness 

il 

Contact ratio 

1< C], C2<2 

Mesh phasing 

pi=0, 0<p2 = h<l 

Variational amplitude 

0< k a i, k a 2 <0.5 


Table 5.1 Parameters of an example system in Figure 5.2. 



Figure 5.3 Instabilities regions when Q, =H 2 =Q,£, — £ 2 — £ ; — analytical solution; 
*** numerical solution. The parameters are from Table 5.1 and c, =c 2 = 1.5, h = 0. 
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The minimum value A ( ^ = 0 is obtained when sin c { 7C = sin c 2 7T = 0 and the unstable region 

vanishes. This is achieved for integer contact ratios cy, c 2 where the number of tooth pairs in 
contact remains constant and mesh stiffnesses are time-invariant. For given, non-integer ci? 

between 1 and 2, A ( ^ is minimized by setting cos(c, — c 2 + 2 h)K — — 1 or c, - c 2 + 2h = ±1, 3 . 

By properly choosing the contact ratios and mesh phasing, the parametric instability regions can 
be dramatically reduced. Figure 5.5 compares the instability regions for three cases. The most 
severe condition for primary instabilities (dash-dot lines) is c, = c, =1.5 and /? = 0, which 

maximizes A^j, in (5.12). This condition is markedly improved by changing the phasing h = 0.5 
so that c, — c 2 + 2/7 = 1 (dashed lines). When the contact ratios are close to integers (solid lines, 
c, = 1 . 1, c 2 = 1.9, h = 0.4 ), the primary instability region becomes even smaller. 

Similar conditions are obtained for secondary and combination instabilities (Lin 
and Parker, 2000c). Unfortunately, the primary, secondary, and combination instability 
regions cannot be minimized at the same time. The conditions reducing the primary 
instability regions (dashed lines, Figure 5.5) result in large combination instability 
regions, and vice versa (dash-dot lines). Depending on specific applications, a trade-off 
may be made to reduce multiple instability regions, though none are true minima (solid 
lines). Adjusting contact ratios and mesh phasing is clearly an effective means to 
minimize instability regions and avoid resonances under operating conditions. 

2. Four-gear Systems: Equal Mesh Stiffness Variations 

Two-stage countershaft systems (Figure 5.1a) have two different mesh frequencies £2j Rf2\ 
which means more instability regions than three-gear systems. We consider the case where the 
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gear facewidth and material are such that the mesh stiffness amplitudes are identical at the two 
meshes (£, =£?=£), although the contact ratios and phasing are not restricted. Depending on the 
ratio R=Z?/Zj, the parametric instability regions associated with Qj and Q may overlap each 
other. For R = m/j (m, j are integers), the l=j instabilities (single mode and combination) of Q\ 
and the I=m instabilities of Q occur simultaneously. Because their instability regions are 
typically the largest, the interactions involving either m or j-1 are of most interest. 



Figure 5.4 Vibration modes for the time-invariant system with parameters in Table 5.1. 
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£ 


Figure 5.5 Comparison of instability regions for various contact ratios and mesh phasing. 
The parameters are in Table 5.1. — c t = c 2 = 1 .5, h = 0.5 ; — c, =1.1, c 2 = 1 .9, h — 0.4; - 
•— c, = c 2 = 1 .5, h-Q. 


When R # m, 1/m for integer m, the 1=1 instabilities from one mesh decouple from 
the l=m instabilities of the other mesh. In this case, instability occurs when IQi or l£2j is 
close to (Op+o^, but these instability boundaries can be calculated independently. For IQ} 
= cq,+ai,+£(J], the condition separating bounded and unbounded solutions is 
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The stability regions associated with each mesh frequency depend on the individual 
contact ratios but are independent of mesh phasing as the two mesh excitations are 
uncoupled. For primary and combination instabilities (1=1), maximum regions occur 
when c,=c,=1.5 and minimum regions require c, =c 2 ={1, 2} . For secondary 
instabilities (1=2), the maximum and minimum conditions are c, =c, ={1.25, 1.75} and 
c, =c 2 ={1, 1.5, 2}, respectively. Figure 5.6a shows the instability regions for R=3/5. 
The primary instability associated with Q 2 — 2(03+ £(72 occurs at Q; = R 21: R(2(t^+ £<7? ) in 
Figure 5.6a. 

When R = m or 1/m for integer m, the parametric excitations from the two meshes 
interact. Consider the case with R=l/m , where the l=m instabilities caused by Qi overlap 
with the primary instabilities caused by Q 2 - Considering instability of the p- th mode 
where Qi=2cOp/m+e<7i and Q 2 =mQi=2(O p +m£(Ji , the boundaries have 

a x =[(D% ) + F"y+(El; ) +G%) 2 ] U2 /(m(O p ) (5.14) 

For example, when £2 2 =20i (R=l/2), the boundaries for secondary instabilities 
(overlapping with 42? primary instabilities) are 

Qt = + O 2 + (^ + GS) 2 ] ,/a (5-15) 

2.CO p 

Figure 5.6b shows instability regions in the (Q\, f) plane for R=l/2. Note the instability 
at Qi~g >3 couples with the instability at and the combined instability region is 

much larger than the case without interaction (Figure 5.6a). Using (5.5) and (5.10) in 
(5.15). the slopes of these boundaries are 
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( 5 . 16 ) 


a \ ~ s\n'(2c^n) + k 2 g2 8 2p sin 2 (c 2 ;r) 

+ 2k gi k g2 S 2 p S 2p sin(2c,;r) sin(c 2 /r) cos(2c, -c 2 + 2 h)x] v2 / 

Minimization of O'/ requires c, = {1, 1.5,2} and c 2 ={ 1,2} for 1 <c i 2 <2. The 

instability region can also be reduced by adjusting the phasing h according to 
cos(2C[ —c 2 + 2 h)K = ±1 with sign the same as sin(2c,;r) . The primary instability regions 
under Qi {1=1) do not coincide with any other instability regions. Other overlap situations 
are possible, such as the Q secondary instability (1=2) overlaps with the Qi fourth 
instability (1=4), but the interaction between these higher instabilities is typically weak 
and the instablity regions are much smaller. Combination instabilities can be analyzed 
similarly. 

3. Three and Four-gear Systems: Unequal Mesh Stiffness Variations 

This general case allows all parameters of the two mesh stiffnesses to differ. In 
contrast with prior cases, the gears may have differing facewidths and material properties 
such that the amplitudes of stiffness variation at each mesh vary independently (Eitt)- 
The contact ratios and mesh phasing are unrestricted. The design of one mesh must 
account for dynamic interactions with the mesh stiffness variation of the other. 

When R * m, 1/m for integer m, there is no interaction between the parametric excitations 
from the two meshes. The Qi instabilities are only affected by £/ and the instabilities 
are only affected by £ 2 . 

When R = m or 1/m for integer m, a mode may be simultaneously driven to 
instability by both mesh excitations. In this case, the first mesh instability regions can be 
significantly affected by the presence of £2 and vice versa. Closed-form boundaries of the 
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form (for primary instability) Q\ = 2<ty+£)<7/+£?<x? for independently varying £ h e 2 are 
cumbersome. Alternatively, simple yet accurate approximations for the instability 
boundaries are obtained by presuming a linear variation of the boundaries in the {Qu £ 1 ) 
plane for given To construct this linear approximation, one point is calculated under 
the condition £ l = 0 ,£ 2 =C and a second point is obtained at e l =e 2 =C . The primary 
stability boundary limits for £, = 0,e 2 = C are 

ft, =2 [(O’+CG™)’] 1 ' 2 (5.17) 

From (5.9), the stability boundary limits for £ x - C,£ 2 = C are 

ft, =2l0 p ±^-[(D“ , r *F";f+(E'»+a^ff 1 (5.18) 

An example is for the primary instability when R=l, C=0.3. Connecting the two 
points obtained from (5.17) and (5.18) yields the instability boundaries, which agree well 
with the numerical solution (Figure 5.7a). Assembling the (i2;, £i ) planes for various £ 2 =C 
generates three-dimensional plots of Qi versus £u £2 (Figure 5.7b). The parametric 
excitation in the second mesh dramatically changes the shapes of the instability regions. 
Notice that the second parametric excitation widens the primary instability region for 
small £\ compared to mono-frequency excitation (Figure 5.3). In contrast, the 
combination instability at £2i~ab +CO 3 disappears near £j=0.23 in Figure 5.7a. In other 
words, an otherwise unstable system is stabilized by the presence of a second parametric 
excitation. The solid line in Figure 5.7b indicates points where the cih +CO 3 combination 
instability vanishes. 
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Figure 5.6 Instabilities regions when Q. x = R£l 2 , £ x -£ 2 = £. (a) R-3/5, (b) R=l/2. The 
parameters are in Table 5.1 and c x =c 2 = 1.5, h = 0. *** denotes numerical solutions. 
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Figure 5.7 Instabilities regions when (a) &i vs. £1 and £?-C-0.3. (b) 12/ vs. £ 7 , £2 

and the solid line indicates vanishing of the combination instability. The parameters are 
in Table 5.1 and c, = c 2 = 1 .5, h- 0 . 
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5.1.3 An Example 

The two-stage gear system (Figure 5.1a) studied by Tordion and Gauvin (1977) and 
Benton and Seireg (1980b) is used as an example. These two papers come to markedly 
different conclusions as discussed below. The system parameters are given in Table 5.1 
and c, =1.47. c, =1.57. In keeping with the published work, the double-tooth contact 
mesh stiffness k max =l is kept constant, and the average mesh stiffnesses k gj decrease as k a , 
is increased (Figure 5.2b). 

Tordion and Gauvin assumed that k v i and k Y j have the same amplitude and 
frequency but different contact ratios and phasing. They applied an infinite determinant 
method (Bolotin, 1964) to plot the boundaries of primary and secondary instabilities 
(dashed lines in Figure 5.8). Their results deviate significantly from the numerical 
solution as a result of analytical errors. In addition, the Fourier expansion they derived for 
rectangular waveforms (equations (11) and (12) in Tordion and Gauvin, 1977) is 
incorrect. Nevertheless, they conclude that “The phase displacement between the 
meshing stiffnesses has a great influence on the width of the instability regions.” 

Benton and Seireg (1980b) considered the same system. They decoupled the 
equations using the modal transformation and neglected the off-diagonal terms of the 
transformed time-varying stiffness matrix (that is, 0 7 K v (/)<J> ). These treatments reduce 

(5.6) to three uncoupled Mathieu equations. The average value of two contact ratios was 
used to make the stiffness variations k v j and k v i identical. With these approximations, 
they conclude that the instability regions are independent of the mesh phasing, that is, 
“the normal mode technique ... without considering the phase variations ... provide(s) a 
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relatively simple means of predicting the instability regions with sufficient accuracy for 
practical purposes.” This conflicts directly with Tordion and Gauvin. In fact, the mode 
uncoupling method does not provide satisfactory results when the mesh phasing is non- 
zero (Figure 5.8b). 

The perturbation results resolve the discrepancy: Mesh phasing strongly impacts the 
mesh stiffness variation instabilities. The excellent agreement of analytical and numerical 
stability boundaries confirms this finding (Figures 5. 3-5. 8). 

To further validate the stability conditions, free responses under non-trivial initial 
conditions are calculated numerically (Figure 5.9) for the parameters at point A of Figure 
5.8 (£2 = 4.2, k a =e=0.3). For point A in Figure 5.8a, the responses are unstable (Figure 
5.9a), as identified by perturbation and numerical methods. This point, however, is stable 
according to Tordion and Gauvin (Figure 5.8a). When the phasing h=0.4 at point A 
(Figure 5.8b), stable responses occur (Figure 5.9b). This is consistent with the 
perturbation and numerical solutions but conflicts with both Tordion and Gauvin s and 
Benton and Seireg's results. 

5.1.4 Discussion 

Rectangular waveforms are close approximations of the mesh stiffness in spur gears 
with involute teeth. For helical gears or spur gears with tooth modification, mesh stiffness 
deviates from the rectangular shape. Equation (5.5) is not valid for other functions, but 
the general Fourier expansion (5.6) can still be used in matrices D, E, F, G to determine 
the instability boundaries. 
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If damping is considered, the system stability improves and the instability regions 
shift to the right in the (£2. e) plane. Furthermore, damping and non-linearity (e.g., tooth 
separation) must be considered to determine the limit cycle amplitude of the dynamic 
response when operating conditions cause instability. 

The instability analysis can be reduced to single mesh gears with one natural 
frequency (ty,. From (5.9), primary and secondary instabilities vanish as the contact ratio 
cj={l,2 } and ci={l, 1.5, 2}, respectively. Maximum primary and secondary instability 
occurs at ci=1.5 and ci={1.25 , 1.75}, respectively. Kahraman and Blankenship (1999) 
experimentally studied a pair of spur gears under mesh stiffness excitation for various 
contact ratios. They showed that the amplitude A/ in the first mesh frequency harmonic of 
the response is minimized when the contact ratio ci={1.0, 2.0}. This is because 
parametric excitations are eliminated for integer contact ratios. When the mesh frequency 
their measured A / reaches maximum for ci~1.4. A possible explanation for the 
high A i at this contact ratio is due to the combined effects of primary and secondary 
instabilities. First, for both primary instability excited by the first harmonic of k(t) and 
secondary instability excited by the second harmonic of k(t), the dominant response 
frequency is a),. We now examine the contact ratios where both instabilities are active. 
Although the maximum primary instability region occurs at C]=1.5, the secondary 
instability region is eliminated there. For cj=1.25, the secondary instability region is 
maximal but the primary instability region is small. For ci~1.4 (average of 1.25 and 1.5 ) 
or 1.6 (average of 1.75 and 1.5), however, both primary and secondary instabilities have 
significant instability regions. Generally, the larger an instability region, the higher 
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response amplitude occurs due to this instability (to see this heuristically, note that both 
the slope of the stability boundaries in (17) and excitation of first order response in (10) 
are proportional to the same quantities D, E, F, and G). Accordingly, for ci~1.4 and 1.6 , 
both instabilities induce large response and jointly contribute to large A /. When the mesh 
frequency f2~2(On (not shown in Kahraman and Blankenship, 1999), the instability is 
caused only by primary instability excited by the first harmonic of k(t). Because the 
primary instability region is maximal at c/=7.5, the response amplitude At around Q~2co„ 
also becomes maximal. Therefore, from the viewpoint of dynamic instability and 
amplitude, contact ratios in the range 1.4-1. 6 are harmful to single-mesh gears at high 
speeds. 

5.2 Planetary Gear Parametric Instability 

Planetary gears have multiple time-varying mesh stiffnesses and their parametric 
instability has not been analytically investigated. In planetary gears, parametric 
excitations are complicated as different contact ratios and phasing conditions exist 
between the sun-planet and the ring-planet meshes. The objective of this study is to 
analyze parametric instability excited by different time-varying mesh stiffiiesses in 
planetary gears. The perturbation method used for two-stage gears is extended to 
planetary gears that have degenerate vibration modes due to their cyclic symmetry. The 
well-defined modal properties of planetary gears are used to derive simple expressions 
for instability boundaries separating the stable and unstable regions. From these 
expressions, the effects of contact ratios and mesh phasing are analytically determined to 
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Figure 5.8 Comparison of instability regions. The parameters are from Table 5.1, 
c,=1.47, c 2 =1.57, and phasing (a) h=0, (b) h=0.4. — Perturbation method; *** 
Numerical method; — Tordion and Gauvin (1977); — Benton and Seireg (1980b). 
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provide insight into planetary gear design. In practice, planet mesh-phasing schemes are 
often applied to cancel or neutralize resonant response at speeds where the mesh 
frequency is near a natural frequency (Seager, 1975; Kahraman, 1994b; Kahraman and 
Blankenship, 1994; Parker, 2000). In this same spirit, this study shows that particular 
parametric instabilities can be eliminated under certain phasing conditions that can be 
achieved by proper selection of design parameters. Tooth separation nonlinearity induced 
by parametric instability is numerically simulated. 


5.2.1 System Model and Mode Properties 

The planetary gear dynamic model used is based on the one developed in Chapter 2 
(Figure 2.2). Translational degrees of freedom in that model are eliminated, and only 
rotational motions of the gear bodies are considered. The sun-planet and ring-planet tooth 
meshes are modeled as linear springs with time-varying stiffnesses k s „(t), kr„(t), n=l, ...A'. 
The system equations of motion are 

Mq + K(/)q = F (5.19) 
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= [u„u r ,u„u l ,-,uj, F = [TJr„ TJr„ TJr„ O.-.Of 


where the summation index n ranges from / to N. In the stiffness matrix K (t), each mesh 
stiffness is represented by k m (t) = k sr + k ln (t), k w (t) = k ip + k 2n (t), n—l,...N , where k sp . 


krp are mean values and ki„, ki„ are time-varying components of the «-th sun-planet and 
ring-planet meshes. For spur gears, rectangular waves are often used to approximate 
mesh stiffnesses alternating between d and d+1 pairs of teeth in contact. Figure 5.10 
shows the mesh stiffness variations ki„, kj„ with peak-to-peak amplitudes 2k vs , 2k vr , 
contact ratios c s , c r , and phasing angles y sn T, (Y sr +Ym)T. The sun-planet and ring-planet 
meshes have the identical mesh frequency Q=2ti/T. The sun-planet mesh stiffnesses 
between planets differ only by a time transition (or phase angle), and likewise for the 
ring-planet meshes. Mesh phasing in planetary gears is determined by planet position 
angles y/„ and the numbers of teeth z s , z r , z p for the sun, ring, and planets (Kahraman and 
Blankenship, 1994). %„= y/^ s /(2 k) denotes the mesh phasing between the first and n- th 
sun-planet meshes ( y sl =0 ); y™ = V^r/(2n) is the mesh phasing between the first and /7-th 
ring-planet meshes (^/=0); y sr =z p /2 is the mesh phasing between the sun-planet and ring- 
planet meshes for each planet. Note for even z p and y sr —l/2 for odd z p . Expansion of 
ki„, fan in Fourier series yields 


k hl (0 = 2k vs £ K? sin IQt + b™ cos /£>/), 
1=1 

k 2 „ (/) = 2* vr £ {a™ sin l£lt + cos lQ.t) 

1=1 


(5.21) 
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a » =— p-sin[//r(c, -2y sn )\sm(ljic s ), a™ =~sin[Itr(c r -2y sr -2y, v )]sin(lnc r ) 

lx lx 

K! = —cos [ln(c s - 2 y m )]sin (/«?, ), ^ = ~-=-cos[//r(c, - 2y lr - 2y„, )]sin(/ac r ) 

in In 

Mesh stiffnesses depend on many parameters including the number of teeth in contact, 

gear facewidth, material properties, profile modifications, and applied load (Kasuba and 

Evans, 1981). Let e x = k vs / k sp and e 2 =k vr /k rp be the relative amplitudes of mesh 

stiffness variation. In the simplest approximation, mesh stiffnesses are assumed 
proportional to the number of tooth pairs in contact, that is, k sp =c s kj h krp=c r ki, and k vs - 
ki,/2 , k vr — ki t /2 , where £/„ ki, are one-pair tooth bending stiffnesses of the sun-planet and 
ring-planet meshes. With this simplifying assumption, the nominal amplitudes of stiffness 
variations are £i=l/(2cj , £2=l/(2c r ), and one obtains the explicit relation £ 2 =£i c s /c r . In 
practice and in the analysis that follows, however, £/ and £2 are not constrained to these 
“nominal” values but vary independently of the contact ratios to account for the many 
factors influencing mesh stiffness variation amplitudes. For this modeling, we let 
£i = £ 2 /g = £ where g=0(l) (and g=c s /c r under the simplifying assumption noted above). 
Substitution of (5.4) and (5.21) into (5.20) yields 

K(/) = K 0 + 2e|](K;, , I ) sin to/ + cos to/) (5.22) 

/=1 

where the average stiffness matrix Kohas the same form as (5.20) with k sn , km substituted 
by k sp , krp. The Fourier coefficient matrices are also in the form (5.20) with k s „, k rn 
substituted by k sp a { ^ , gk ^ for K^and substituted by k sp b ( " , gk ^ for K ( v ' 2 \ 

For the time-invariant case, the eigenvalue problem associated with (5.19) is 
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K 00, = <y, 2 M 4> I 


(5.23) 


The natural frequencies and vibration inodes have unique properties (3.7), (3.14), (3.20). 

We specialize these properties for the case of a rotational vibration model with fixed ring 

and L=N+2 degrees of freedom. In this case, all vibration modes can be classified into 

one of three categories: (1) a rigid body mode ( O)i=0 ), (2) two modes with distinct natural 

frequencies ( cih , (Ql), and (3) a group of degenerate modes with multiplicity N-l 

(ft£=...=tf>v /). In the distinct modes, all planets have identical motion 

u„= ui, n=l,2 N (5.24) 

In the degenerate modes, the carrier, ring, and sun have no motion, and the planet 

rotations satisfy 

u„ = [M,sin(^ 2 -^„) + w, sin^J/sin^, n=l,2 N (5.25) 

These well-defined properties are valid not only for equally spaced planets with position 
angles y/„ =27t(n-l)/N, but also for diametrically opposed planets with y/„+N ? =y/„+n. 

5.2.2General Expression for Instability Boundaries 

We determine the operating conditions (that is, mesh frequency Q and stiffness variation 
amplitude e) that lead to instability when IQ-cop+a), for integer /. The rigid body mode 
( 0 ) 1 = 0 ) is not excited under operating conditions, and does not affect the instabilities of 
other modes. Only the two distinct modes and the group of degenerate modes are 
considered in what follows. The parametric instability when IQ is close to (Op+a \ is 
considered. Let IQ = a> p + (O q + £<T , where <7 is a detuning parameter to be determined. 

From (5.8), when (tip and (fy are both distinct, the instability boundaries are 


N AS A/CR— 200 1-210939 


85 





Figure 5.10 Modeling of mesh stiffness variations in the sun-planet and ring-planet 
meshes. c s , c r are contact ratios, and y s „, y™ , Ysr are mesh phasing. 


single mode instability: £2 = 
combination instability: £1 = 



(5.26) 

(5.27) 


where A ( ^ = {D^j 2 + (E ( ^) 2 . 

When (Op or o\, are degenerate, the complexity of instability solutions depends on 
the multiplicity of the degenerate natural frequencies. We first study the case with 
multiplicity two, say (Op=a The single-mode instability boundaries for (Of,=OJp are (Lin 
and Parker, 2000d) 
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When (tyrcty are degenerate but o\ is distinct, the combination instability is 


£2 = ± £[(Z)«D«' + + £«£"’ + < 5 - 29 > 

When the degenerate natural frequencies =C 0 m *2 have multiplicity m > 2, the 

critical <7 is obtained by requiring the real parts of the eigenvalues of the following matrix 


to be non-positive, 


e£ d«-(o»,)i 


3 ntxm 


I>i!i +(<76> 3 )Imx BI 


-E 


(O 


(5.30) 


where D^.E^ are mxm submatrices of D , E' including rows and columns 3, ..., 
m+2. Generally, no closed-form solution can be derived for these single-mode 
instabilities, but numerical evaluation can determine <7 and the instability boundaries 
unless have special features (e.g., diagonal matrices) as discussed later. More 

is possible for combination instabilities. For combination instability of distinct (O and 
degenerate (O ^- ... =&V 2 , the instability boundaries are 

(5.31) 


a = ±7 + CO'*®-®.) 

* * <7=3 


The above expressions apply for a general system with degenerate natural 
frequencies. These results reduce to simple forms when specialized to planetary gears. 
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5.2.3. Planetary Gear Parametric Instability 

Planetary gears’ well-defined modal properties are now used to simplify the above 
instability conditions to compact expressions suitable for use in applications. Mesh 
stiffness are approximated by rectangular waveforms (5.21). Expansion of D (, \ E <]> gives 

N N 


D (,) =(p T iC'ty =k Y a 0) 8 {p) d { 

pi / *p vl rq sp sn s\ 


q) 



n-\ w=l 

N N 



(5.32) 


n= 1 


n = 1 


where S[ p) = - u c + u n is the deformation of the «-th sim-planet mesh in mode and 


8 ( n p] ] = —u c - u„ is the deformation of the «-th ring-planet mesh. For different spacing and 

phasing of planets, matrices D ffl , E (1> have special features that simplify the instability 
conditions. 


1. Equally spaced planets 

For equally spaced planets (y/ n =2n(n-l)/N), the vibration modes have structured 
properties (5.24) and (5.25). In this case, (z s +: r )/N = integer and the planet meshes have 
phasing y sn =(n-l)z s /N and y r „=(n-l)z r /N, n=l,...,N. The two possible phasing conditions 
are examined below. 


(1) In-phase Planet Meshes 

Consider the case when all the sun-planet meshes are in-phase ( y s „=0 , n= 
and all the ring-planet meshes are in-phase (y m =0, n=l,...,N). This design is typical when 
optimizing load sharing and results when z s and z r are each integer multiples of N. 
However, there is a constant phasing y sr =0 (even z p ) or 1/2 (odd z p ) between the sun- 
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planet and ring-planet meshes for each planet. The Fourier coefficients o { ’^ ,b['J MJ 

in (5.21) are independent of the planet index n. 

For a distinct natural frequency vibration mode property (5.24) leads to 
= 8^f } , 8^ = 8j:( >) for any planet n. Using (5.21) and (5.32) on (5.26), the primary 


instability boundaries have 




Tt 


(5.33) 


2 k sp k rp g(5 ( s i ) 8 ( r i ) f sin(^c 5 )sin(/ar r )cos^'(c J —c r + 2y sr )] 


Obviously, A" = 0 when c s , c, are integers and all instabilities vanish. If c s , c r * integer, a 

second choice to reduce the instability regions is to set the third term in (5.33) to be 
negative by adjusting c s , c r , and y sr . From (5.27), the combination instability boundaries 
for two distinct modes have 


A ( ^ = (—) 2 [(K p 8 ( sl p) 8 { f sin 7cc s ) 2 +(k n ,g8 ( f8lf sin nc r f + 


n 


(5.34) 


2 K p k, v g8^8f8\?8f sin(^ s )sin(^ r )cos^(c s -c r + 2 y sr )] 


Inertias (kg) 

I c /r 2 =6, l s /r 2 =2.5, I p /r p 2 =2 

Planet mass (kg) 

> 

II 

Mesh stiffness (N/m) 

II 

II 

Natural frequencies 
(kHz) 

a)i=0, ob=1.212, 
0)3=0)4=1.592, 0)5=2.196 


Table 5.2 Parameters and natural frequencies of an example planetary gear with fixed 
ring and three planets. 
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For the degenerate natural frequencies 003=. ..=004,1, the single-mode instability 
boundaries are the same as (5.28) with 

r m =(£%)’ + <C > 2 =(-^) ! [(*,sin/®,) ! +(k T g S mlnc,y- + 

2k sp k ip gsm(hte s )sin(lm' r )cosln(c s -c,.+2y sr )] 

Equation (5.35) applies to any number of planets N. 

For the combination instability of degenerate (05=. ..=<%+/ and distinct oy, the 
combination instabilities (/=/) of a distinct mode and a degenerate mode always vanish 
for any N when the planet meshes are in-phase. A more general condition for vanishing 
of these combination instabilities is that each of (lz s ±l)/N, (h r ±l)/N * integer. 

As an example, Figure 5.11 shows the instability boundaries for a planetary gear 
with three equally spaced planets. The parameters and nominal natural frequencies are 
given in Table 5.2 and the vibration modes are shown in Figure 5.12. The mesh phasing 
is Ysn~Ym~0, y sr =]/2 and the contact ratios are c s =1.4, c r =1.6. We specify £j=£ 2 =£ (that 
is, g=l). The natural frequencies change as ovaries because the average mesh stiffnesses 
depend on the amplitude of the mesh stiffness variation. The analytical solutions (solid 
lines) from (5.26), (5.27), (5.28), and (5.29) agree well with the numerical solutions using 
Floquet theory and numerical integration. Note the combination instabilities at 002 + 0 ) 3 , 
( 03 + 0)3 vanish because they involve the distinct (00,005) and degenerate ( 003 = 0 ) 4 ) natural 
frequencies. Figure 5.13a shows the primary instability regions for different contact 
ratios c$, cy; the stiffness variation amplitude £=0.3 and phasing ^ r =0. All instabilities 
vanish when the contact ratios are integers. The size of the 200 instability region is 
primarily affected by c r while insensitive to changes in c s . This is because the dominant 
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deformation of mode <jh occurs in the ring-planet meshes (Figure 5.12). The 2(0 . 5 
instability region is mostly affected by c s because most deformation occurs in the sun- 
planet meshes for this mode. The primary instability of the degenerate modes at 2 0)3=20)4 
is influenced by both c s and c r . The primary instability regions are maximized for 
c.f=c r = 1.5, which is expected from (5.33) and (5.35). The relative phasing y sr between the 
sun-planet and ring-planet meshes can have a major impact. When the phasing y sr =l/2 
and other parameters are the same as in Figure 5.13a, the 20)3=20)4 instability vanishes 
for any c s =c r (Figure 5.13b). This is because T f/) =0 in (5.35) for k sp =k rp , g=l and c s =c r . 

The above analyses show that contact ratios and mesh phasing significantly affect 
the operating condition instability regions. In practice, particular instabilities can be 
minimized by proper selection of contact ratios and mesh phasing, which are adjusted by 
center distance, diametral pitch, pressure angle, tooth addendum, numbers of teeth, and 
other parameters. 


(2) Sequentially Phased Planet Meshes 

Here we consider equally spaced planet systems where the sun-planet and ring- 
planet meshes are sequentially phased with %„=(n-l)z s /N and y rt ,=(n-l)z r /N. This case 
corresponds to r, /N, z r //V> integer but (z s +z r )/N= integer. A constant phasing y sr exists 
between the sun-planet and ring-planet meshes for each planet. For the sequential 
phasing, the Fourier coefficients in (5.21) satisfy 


Y N a 0) = T N a U) = y N b 0) = y N b 0) = 0 


(5.36) 


when lz s /N, lz r /N are non-integer. 
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For distinct natural frequencies (q,. the primary instabilities (1=1) of distinct modes 
always vanish when the planets are sequentially phased. Physically, it means that these 
instabilities are not excited because the resultant modal excitations from the sun-planet 
and ring-planet meshes each are zero. The secondary instabilities (1=2) of distinct modes 
also vanish when 2z s /N , 2: r /N are non-integer. Similarly, combination instabilities (1=1) 
of distinct ay and ay, D ( ^ = E™ = 0 always vanish in this case. 

For instabilities involving the degenerate modes, it is difficult to obtain simple 
expressions for the instability boundaries; their instability conditions can be calculated 
from (5.28), (5.30), and (5.31). A special case is the combination instability (1=1) of 
degenerate and distinct ay when (z s ±l)/N, (z r ±l)/N each is non-integer. In 

this case, the combination instabilities of distinct and degenerate modes vanish. 

Figure 5.14 shows the instability boundaries for the same system as in Figure 5.1 1, 
except the three planets are sequentially phased with y s „=[0, 1/S, 2/3], y m =[0, 2/3, 1/3]. 
Note the primary, secondary, and combination instability regions vanish for the distinct 
natural frequencies ay, ay; only instabilities involving at least one of the degenerate 
modes ay=ai) exist. 

The foregoing results for equally spaced planets are summarized in Table 5.3. 
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Figure 5.1 1 Instability regions for the system in Table 5.2 and in-phase meshes y„-ym 
y sr =l/2. c s =1.4, c r =1.6, £=£i =£2 (g=l)- — analytical solution; *** numerical solution. 



o> 2 =121 2 Hz 



o> 3 =1592 Hz 




co 4 =1592 Hz 


Figure 5. 12 Mode shapes of the system in Table 5.2. The carrier motion is not shown. 
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Figure 5.14 Instability regions for the three-planet system in Table 5.2 and sequentially 
phased planet meshes with y s „=[0, 1/3, 2/3], ^=[0, 2/3, 1/3]. c s =1.4, c r =1.6, £=£;=£: 
(g=l), y sr =l/2. — analytical solution; *** numerical solution. 


Planet Mesh 
Phasing 

Single-Mode Instabilities 

Combination Instabilities 

Distinct Mode 

Degenerate 

Mode 

Distinct + 
Distinct Mode 

Degenerate + 
Distinct Mode 

In-phase 

— = integer 

N N 

from (5.26) 

from (5.28), 
(5.30) 

from (5.27) 

always vanish 

Sequentially phased 
— * integer 

N N 

primary always vanish; 
secondary' from (5.26) 
and vanish if 2z s /N, 
2z r /N ± integer 

from (5.28), 
(5.30) 

always vanish 

from (5.31) and 
vanish if (z s ±l)/N, 
(Zr±l)/N # integer 


Table 5.3 Instability boundary solutions when the planets are equally spaced (satisfying 
(: s +z r )/N = integer). 
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2. Unequally spaced planets 

When the planets are arbitrarily spaced, the structured modal properties do not exist 
and the general expressions for instability boundaries cannot be further simplified. For 
the practically important case of diametrically opposed planets, however, the vibration 
modes retain the well-defined properties (5.24) and (5.25). The following discussion 
focuses on this case for N/2 pairs of diametrically opposed planets. For the sun-planet 
meshes, each pair of diametrically opposed planets are in-phase (%(„+N 2) = Y») for even z s 
and are counter-phased (y S („+N2) = Ysn + l/Z) for odd z s . Analogous rule applies for the ring- 
planet mesh phasing. Note that adjacent planets have arbitrary mesh phasing Ysn = Wr^r 
/(2n), y r „=y/^ r /(2K), n=l,...,N/2. 

When the sun-planet and ring-planet meshes are both counter-phased (odd z s , z r \ 
equation (5.36) holds for odd / because a { £ = -a^ n+NI2) and similar relations for 

°ra’ Kn- Recalling modal property (5.24), =0 in (5.32) for distinct (Op, 

(fy and odd /. From (5.26) and (5.27), primary and combination (/=/) instabilities of 
distinct modes always vanish in the counter-phased case. Physically, these instabilities 
are eliminated because the modal excitations from each pair of diametrically opposed 
planets always cancel each other. The instability regions involving degenerate modes are 
obtained from numerical evaluation of the eigenvalues of (5.30). 

When pairs of opposing sun-planet or ring-planet meshes are in-phase (even z s or 
even : r ), no simple expressions for instabilities regions are available; numerical solutions 
are obtained from (5.26), (5.27), (5.30), and (5.31). 
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5.2.4 Dynamic Response and Contact Loss 

When planetary gears are operated inside an instability region, damping and 
nonlinearities from friction, tooth separation, etc bound the unstable linear model motion. 
Figure 5.15a shows the RMS steady-state planet response amplitude versus mesh 
frequency for the same system as in Figure 5.11 and the stiffness variation £=0.3. 
Rayleigh damping C=(0.07)*M+(0.07)*K is added to system (5.2) and the force vector 
is F =[-2000 1000 0 0 0] T N. The solutions are obtained from numerical integration using 
mesh stiffnesses in rectangular waveforms (Figure 5.15b). The degenerate modes ( 03 = 0)4 
only have planet motion, so they are not excited for the torques applied to the carrier and 
sun. The combination instability 00 + 0)5 is sufficiently damped so that it is not apparent 
in Figure 5.15a and 16a. The resonance excited by the primary instability 2a*, is 
extremely large because tooth separation is not considered; the mesh stiffnesses are pre- 
specified functions of time (Figure 5.15b). In practice, tooth separation (clearance 
nonlinearity) occurs for large dynamic responses and its effects are dramatic. Figure 
5.16a shows the response for the same system as in Figure 5.15a, but tooth separations is 
modeled. The mesh stiffness k s „ or k r „ is set to zero if the corresponding tooth 
deformation S s „<0 or 4,<0 at any step of the integration. The response amplitude of the 
2o)5 primary instability is significantly reduced from that in Figure 5.15a. Moreover, a 
softening jump phenomenon occurs. Sun-planet tooth separation ( k sp =0 ) is apparent in 
Figure 5.16b for £ 2 = 2 ( 05 . The interactions of mesh stiffness variation and clearance 
nonlinearity was studied by Kahraman and Blankenship (1996, 1997) for single-mesh 
gears. Their effects on multi-mesh planetary gears need additional investigation. 
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Figure 5.15 (a) Steady state RMS of dynamic planet rotation versus mesh frequency 
Q when tooth separation is not considered. The parameters are as in Figure 5.11 with 
£= 0.3 . (b) The sun-planet mesh stiffness is pre-specified as shown. 
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Figure 5.16 (a) Steady state RMS of dynamic planet rotation versus mesh frequency fJ 
when tooth separation is considered. The parameters are as in Figure 5.1 1 with £=0.3. 
Circles (o) indicate increasing speed Q and crosses (x) indicate decreasing Q. (b) Sun- 
planet tooth separation ( k sp =0 ) occurs for Q=4.5kHz=2(Os. 
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CHAPTER 6 


SUMMARY AND FUTURE WORK 


6.1 Summary and Benefits 

This work analytically investigates several key issues in planetary gear dynamics. 
The main results and their benefits are summarized for each specific topic. 

1. Dynamic Model of Planetary Gears 

A lumped-parameter model is developed for spur planetary gears. Critical factors to 
gear vibration are considered, including arbitrary planet spacing and phasing, gyroscopic 
effects that are important in high-speed applications, mesh stiffness variation, and 
transmission error excitation. The model is applicable to general epicyclic gears with 
various configurations. It is suitable for the dynamic analysis of critical issues in 
planetary gear vibration and useful for design guidance (response, natural frequencies, 
etc.). This model is a key building block that can be expanded to couple with the housing 
and include multiple stages. 

2. Modal Properties 


NAS A/CR— 200 1-210939 


101 



This study characterizes the natural frequency spectrum and vibration mode 
properties in planetary gears. Planetary gears possess rich modal structure that is crucial 
for eigensensitivity, dynamic response, and stability analysis. Rigorous characterization 
of these special properties is a fundamental advance for planetary gear vibration research. 
These properties are capsulized as follows: 

• When N planets are equally spaced, the vibration modes can be classified into six 
rotational modes with distinct natural frequencies, six pairs of translational modes 
with degenerate natural frequencies of multiplicity two, and three groups of planet 
modes with degenerate natural frequencies of multiplicity N-3. Each type of vibration 
mode has unique properties due to the cyclic symmetry of the system. 

• When N/2 pairs of planets are diametrically opposed, the rotational and planet modes 
have the same structure as for equally spaced planets; translational modes lose their 
degeneracy but retain their distinct properties. 

• When the planets are arbitrarily spaced, the rotational and translational modes 
generally lose their well-defined structure; the remarkable properties of planet modes 
are not affected by planet spacing. 

3. Design Parameter Variations 

The effects of design parameter variations on planetary gear free vibration are 
analytically investigated. The sensitivities of natural frequencies and vibration modes to 
key parameters are expressed in simple, closed-form formulae that allow one to quickly 
identify the parameters that most impact the modal properties. Well-defined veering rules 
are derived to examine the dramatic changes of vibration modes that are possible even for 
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small parameter variations. These results rely on the structured vibration mode properties 
characterized before. For design use, the results allow qualitative and quantitative 
assessment of the effects of design parameters on planetary gear free vibration: 

• Each of the sun-planet and ring-planet mesh stiffnesses affect one rotational mode, 
one translational mode, and one planet mode. 

• The transverse support stiffnesses and component masses of the carrier, ring, and sun 
only affect translational modes. 

• The torsional support stiffnesses and moments of inertia of the carrier, ring, and sun 
only affect rotational modes. 

• The planet bearing stiffness and inertia are critical parameters that affect most modes. 

• The operating speed does not affect rotational and planet modes; translational mode 
natural frequencies at zero speed split into distinct ones as speed increases. 

4. Mesh Stiffness Variation Excitation 

The parametric instabilities excited by mesh stiffness variations are investigated for 
two-stage and planetary gear trains. The operating conditions leading to parametric 
instabilities are analytically determined and numerically verified. The following findings 
are obtained from this study. 

• The contact ratios and mesh phasing strongly impact the instabilities induced by mesh 
stiffness variation. Simple, exact formulae derived from perturbation analyses provide 
design guidance to suppress particular instabilities by adjusting the contact ratios and 
mesh phasing. 
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• For two-stage gears in a countershaft configuration, the mesh excitations interact 
when one mesh frequency is an integer multiple of the other and dramatically change 
the instability conditions compared to two decoupled excitations. 

• For planetary gears, the structured modal properties lead to vanishing of certain 
instabilities, including the combination instability of distinct and degenerate natural 
frequencies when planet meshes are in-phase, and the single-mode and combination 
instabilities of distinct natural frequencies when planet meshes are sequentially 
phased. Along with contact ratios, selection of planet phasing is an effective means to 
avoid parametric instability and its associated large response. 

6.2 Future Work 

To thoroughly understand planetary gear dynamics and develop reliable design 
tools for noise and vibration reduction, this research needs to be advanced in several 
challenging areas discussed below. 

1. Nonlinear Effects 

Tooth separation nonlinearity strongly affects the dynamic behaviors of single-mesh 
gears, including jump phenomena and multiple steady state solutions (Blankenship and 
Kahraman 1995; Kahraman, 1992; Kahraman and Blankenship, 1996, 1997; Kahraman 
and Singh, 1991; Parker et al., 2000b; Rook and Singh, 1995; Theodossiades and 
Natsiavas, 2000). Typical methodologies such as direct integration, harmonic balance, 
perturbation, and shooting method were used in the above analyses. This nonlinearity has 
not been fully investigated in planetary gears. The contact loss nonlinearity, coupled with 
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mesh stiffness variation, transmission error excitation, and mesh phasing, is expected to 
impact dynamic responses, load sharing among planets, bearing loads, tooth fatigue, and 
gear noise. Analytical study of nonlinear effects will provide improved design guidance 
for planetary gears based on dynamic response. 

2. Transmission Error Modeling 

For single-mesh gears, static transmission error (STE) is naturally defined and 
widely used as a metric for vibration and noise (Gregory et al., 1963-64; Smith, 1987; 
Ozguven and Houser, 1988b). The specification of transmission errors in planetary gears 
is complicated due to the carrier rotation and multiple meshes coupled together. Donley 
and Steyer (1992) simply used a “net transmission error” that is the average of all 
individual transmission errors of the sun-planet and ring-planet meshes. Although this 
simplification reduces the computational and modeling effort, it is a significant 
simplification difficult to interpret of the dynamic mesh forces. Other studies (Kahraman, 
1994b, 1999; Kahraman and Blankenship, 1994) used Fourier series to represent STE in 
each mesh, but did not address the physical modeling related to profile error, pitch error, 
run-out error, and misalignment. It remains a question how to synthesize transmission 
errors from these parameters in planetary gears. More investigations are needed to define, 
calculate, and measure a useful analog of STE in planetary gears. 

3. Dynamic Response 

Dynamic response from analytical models is useful to predict the noise and 
vibration early in the design stage. To reliably determine dynamic response, it requires 
precise representation of dynamic excitation from tooth meshes. The time-varying mesh 
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stiffness and transmission error for single-mesh gears can be calculated from the contact 
analysis program LDP (Houser, 1990) or other tools. A major challenge is to match 
dynamic response obtained from analytical models using these approximations as input 
with benchmark computational and experimental results. No such comparison exists for 
planetary gears, and this restricts the application of analytical models in practical design. 
In addition, the effects of major design factors ( e.g ., contact ratio, mesh phasing, mesh 
stiffness variation, support/bearing stiffness, tooth modification) on dynamic response 
need be characterized for design guidance. This investigation will start with two-stage 
spur gears and then extend to more complicated planetary gears. 

4. Ring Gear Flexibility 

The ring body is considered rigid in this dissertation. In practice, internal gears with 
thin rims may distort elastically. Experimental measurements (Ma and Botman, 1984) 
and finite element computations (Kahraman and Vijayakar, 2000) have shown dramatic 
changes in load sharing among planets from the ring flexibility. Because the ring is 
usually connected to the housing, its flexibility directly influences the transmission of 
structure-borne noise. The elastic distortion of the ring can also affect the mesh 
stiffnesses and lead to contact loss. Qualitative and quantitative analyses of the effects of 
ring flexibility is necessary to improve understanding of physical behaviors. The 
fundamental task is to extend the current discrete model by including the continuous 
deformation of the ring body. The analytical and computational methods discussed in 
this work can likely be used on this extended model. 

5. Experimental Verification 
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5. Experimental Verification 

The inherent complexity of gear dynamics requires experimental validation of 
analytical and numerical results. For single-mesh gear models, there are extensive 
experiments verifying analytical findings as well as presenting questions for research and 
design. The analytical models of planetary gears, however, have not been properly 
validated because of the lack of benchmark experimentation. The calculated vibrations 
from the finite element tool Calyx agree well with measured results at NASA Glenn 
Research Center (Krantz, 1992) and General Motors (Blankenship and Kahraman, 1996; 
Kahraman and Blankenship, 1997, 1999a, b; Parker et al„ 2000b). More experimental 
studies are needed to develop reliable dynamic models suitable for practical design and 
investigate critical issues in planetary gear vibration. 
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APPENDIX A: SYSTEM MATRICES 


M,, = diag(m h . m h , / r A 2 ), h = c,r,sX---,N 


G = 


0 - 2 m h 0 


2m h 0 0 


0 0 


K n = diag(m c , m c , 0 , m r ,m r ,0, m s , m s ,0,m p ,m p ,0,---, m p ,m p ,0 ) 


K 6 =^(K fA ,K rft ,K, ft ,0,-,0) 


II 

<1 

* 

diag(k hx A, n ,k hll ) 

1, 

h = c, 

r,s 



Z K "i 0 

0 

<2 

k 2 c2 - 





0 

k‘ 2 

K; 2 • 

TS N 

IV 


K m = 

Ik 

n i/ 1 
jl ^52 

k; 2 

0 

•** K^ 2 

0 



symmetric 







- 




X 



10 - sin y n 




[ 

II 

e 

1 cos 





\C n - 
K c2 “ 


symmetric 1 




L 


K" 3 = diag(k ptl k pn 0) 
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sin y/ n 0 

- cos y/ n 0 
-1 0 



sm ~ Vrn - cos yr n sin y/„, - sin y/ n , 

•C, =k sn (t) cos V™ cos^„, 

symmetric 1 

- sin if/ ni sin a r - sin \f/ rn cos a r sin y/ rn 
K r 2 =^r«( / ) cos^sina,. -cos^ n7 cosa r -cosy/ rn 
sin a r -cosa r 1 

sin 2 a; - cos sin cr,. - sin 

K” 3 = £ r „(0 cos 2 or r cosa ; . 

symmetric 1 

sin 2 y/ rn - cos \f/ sn sin y/ sn - sin yr m 

K"i = *„(0 COsV,„ COS ^vn 

symmetric 1 

sin ^ n sinor s sin y/ s „ cos a s - sin y/ sn 

k ”2 = *«(0 -cos^sino', -cos^ Jn cosa s cos^„ 

- sina s - cosa s 1 

sin 2 a s cos a s sin a s - sin a s 

K .o = *m(0 COS 2 or, -cos a, 

symmetric 1 

T(o = [oa t c ,0,0, r r ,o,o, 7; a • • • of 

F(O = [0,F r ,F J .F,-FJ r , 
r = * ra e ra [sin^ m ,-cos^ m ,l] 7 , 
r = k sn e sn [sin Vsn ~ COS \ff „ ,1]^ 

F n =^ m [s in «r + * 4 B e J „[sma J ,-sina r ,l] r 

Note that the mesh stiffness variation is modeled by the time-varying stiffnesses 
k s „(t) and k m (t). 
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APPENDIX B: EXAMPLE SYSTEM 


The planetary gear used in the U.S. Army’s Helicopter OH-58 is the benchmark 
example in the work. Two cases of the system are used in the study. 

Case I: The ring gear is free and the system has 3(N+3) degree of freedom, where N is 
the number of planets. The model parameters are given in Table C.l. 

Case II: The ring gear is fixed ( k r = k ru = 0) and other parameters are the same as Table 
C.l. The system has four planets and 18 degrees of freedom. 



Sun 

Ring 

Carrier 

Planet 

Mass (kg) 

0.4 

2.35 

5.43 

0.66 

I/r 2 (kg) 

0.39 

3.00 

6.29 

0.61 

Base diameter (mm) 

77.42 

275.03 

177.8 

100.35 

Teeth number 

27 

99 


35 

Mesh stiffness (N/m) 

Kp = krp= k m = 5xl<f 

Bearing stiffness (N/m) 

II 

Cl 

II 

It 

II 

Torsional stiffness (N/m) 

k ni = I (f k su = k cu = 0 

Pressure angle (°) 

0Cs = a, = a = 24.6 


Table B.l Model parameters of the planetary gear in the U.S. Army’s helicopter OH-58 
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